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Abstract 

One may care to study operators and their inverses under conditions for which poly- 
nomial decay is too slow, and exponential decay is too fast or not exactly preserved. 
In this paper I establish that a broad class of Banach *-algebras of infinite integral 
operators, defined by the property that the kernels of the elements of the algebras 
possess subexponential off-diagonal decay, is inverse closed in i3(L^(]R'^)). This means 
that each subalgebra of i3(L^(]R'^)) under consideration contains the inverse of each 
of its invertible elements. 

A dividend of the above work concerning symmetry of Banach algebras is demon- 
strated en route to proving the above facts about decay of kernels of integral operator 
inverses. A *-algebra is symmetric if the (algebraic) spectrum of positive elements is 
positive. Historically it has often been quite difficult to verify that a given involutive 
algebra is symmetric, but using techniques discovered by Grochenig and Leinert I will 
show that the algebras under consideration are symmetric. 

In the second part of this dissertation, I present the results of the IEEE Transactions 
on Communications paper written jointly by Thomas Strohmer and me. We develop a 
comprehensive framework for the design of orthogonal frequency-division multiplexing 
(OFDM) systems, using techniques from Gabor frame theory, sphere-packing theory, 
representation theory, and the Heisenberg group. We begin with a Gaussian pulse, 
which has optimal time-frequency localization (TFL), and orthogonalize it (using 
a clever method first discovered by the chemist Lowdin) with respect to the time- 
frequency shifts of an OFDM system. This process retains the optimal amount of the 
TFL of the Gaussian, vastly improving the robustness of the standard OFDM system 
to dispersion caused by multipath propagation and the Doppler effect. Our next 
contribution to the OFDM algorithm is the introduction of Lattice-OFDM, wherein 
we admit more general lattices than the usual separable lattices. This provides even 
more defense against dispersion, and we present a theorem giving the OFDM system 
designer the ability to optimally adapt the lattice and pulse to the channel conditions. 
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Chapter 1 



Introduction 



1.1 Banach algebra results 

This dissertation is a study of off-diagonal decay properties of kernels of invertible 
integral operators and their inverses. Previous strides in the area have been made 
by Gel'fand, Grochenig, Leinert, Jaffard, and Gohberg, among others (see Section 
1.1.2). I extend the results in [23] on the off-diagonal decay of infinite matrices and 
their inverses, and the symmetry of infinite matrix algebras, to the appropriate class 
of distributions on M^*^. 





Figure 1.1: Kernels, of integral operators, exhibiting off-diagonal decay 



One is often compelled to study operators and their inverses under conditions for 
which polynomial decay is too slow, and exponential decay is too fast or not exactly 
preserved [23] . My result establishes that a broad class of involutive Banach algebras 
(Banach *-algebras) of infinite integral operators, defined by the property that the 
kernels of the algebra elements possess subexponential off-diagonal decay, is inverse 



closed in i3(L^(M'^)). This means that each subalgebra of i3(L^(]R'^)) under consider- 
ation contains the inverse of each of its elements. 

An additional result, concerning symmetry of Banach algebras, is established en route 
while proving the above facts about decay of kernels of integral operator inverses. 
An involutive algebra is symmetric if the spectrum of positive elements is positive. 
Historically it has often been quite difficult to verify that a given involutive algebra 
is symmetric. Losert recently solved a 30-year-old conjecture that the group algebra 
of a compactly generated, locally compact group of polynomial growth is symmetric. 
The result of Grochenig and Leinert includes the establishment of symmetry for a 
large class of bi-infinite matrix algebras. Extending their result about the inverse 
closedness of Banach algebras of bi-infinite matrices to inverse closedness of unital 
Banach algebras of infinite integral operators permits me to conclude that the algebras 
I consider are symmetric as well. 

1.1.1 Applications: 

1. Continuous-time blurring by convolution. It is recognized that implementation 
will require discretization, quantization, and truncation, but a thorough understand- 
ing of the continuous theoretical model is often crucial to avoid setting up the model 
in a non-optimal way, of which there are examples from non-uniform sampling theory. 

2. Wireless channel scattering functions. In practical terms a channel scattering 
function is often the sum of a delta "function" and a compact operator^, or the sum 
of a delta and a low-rank operator, each of which is invertible. 

3. Any integral equation whose kernel has appropriate decay properties but is not 
singular. 

1.1.2 Previous results in the field 

To illustrate the context of my main result, I briefiy discuss previous results in the 
area, all of which regard operators on £^(Z'^): 

1. Commutative Cases: 

(a) Wiener's Lemma (Wiener, illustrated in e.g., [44]) : i^{Z'^) is an algebra in which 
the multiplication is convolution. For a G £^(Z'^) consider the convolution operator 
Ta acting on £^(Z^), the action of which is given by TaC = a * c. Wiener's Lemma in 
this case is: if a G i^{'Ij'^) and if T^ is invertible as an operator on £^(Z'^), then a is 



^This is the case when the channel is described by Ricean fading, where there is a line-of-sight component in the 
received signal. 
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invertible in £^(Z'^) and hence the inverse operator is of the form T~^ = T^ for some 
b G e{Z'^). 

(b) Weighted Wiener's Lemma (imphcitly Gel'fand, Raikov, Shilov [20], exphcitly 
Strohmer [53]): Let A = {an} be a Hermitian positive definite Laurent matrix, 
generated by the 0*^ row denoted ak, with inverse A^'^ = {bke}, generated by the O*'' 
row denoted bk, and let v{k) be a positive submuhiphcative {v{k + i) < v{k)v{tj) 
weight function such that Xlfegz \^k\'^{^) < ^^- If "^ satisfies the condition 



hm — -^=^= = 1 and hm yvl 

n-.oo y^(_^) n-*oo 



then Y.kez\M^W < '^■ 

2. Non-commutative cases: 

(c) Non- commutative Wiener's Lemma (Sjostrand [51]): Let C = Ck/ , k,£ & Z'^, 
be a matrix with \Ck/\ < a{k — i) for some a G £^(Z'*), i.e. C is dominated by 
the Laurent function generated by a (see Section 2.6). If C has a bounded inverse D 
on £^(Z'^), then there exists b G £^(Z'^) depending only on a such that \Dk,e\ < b{k — £). 

(d) Exponential off-diagonal decay of infinite matrices (Jaffard [26], Gohberg [21]): 
Suppose A is bounded on i (Z.'^), with bounded inverse. If either ake = for 
\k-i\> Morif lofc^l = C(e-°l'=-^l)> " > 0, then 5 = A^^ satisfies la^l = C(e-^l''-^l) 
for some j3 > 0. 

(e) Polynomial off-diagonal decay of infinite matrices (Jaffard [26]): li A is boundedly 
invertible on £^(Z'^) and \ake\ = 0{\k — ^|^^) for some s > d, then B = A^'^ has the 
same polynomial-type off-diagonal decay \bke\ = 0{\k — i\~^). 

(f) Sub- exponential off-diagonal decay of infinite matrices (Grochenig and Leinert 
[23]): This case includes (c). Let denote the Hadamard product and let Vki : = 
v{\k — i\). li A is boundedly invertible on £^(Z'^) and the row- and column-sums of 
\A v\ are bounded, then B = A~^ has the same off-diagonal decay as A (i.e., the 
row- and column-sums oi \B Q v\ are bounded). Or, if \ake\ < C v~^{\k — i\) for some 
C, then \bke\ < C'v-^{\k - e\). 

The contribution to the field provided by the first part of this dissertation is the 
extension of the results of Grochenig and Leinert to the continuous case. That is to 
say, I will show that if K is an integral operator on L^(]R'^), and the kernel k oi K 
has sub-exponential off-diagonal decay measured either pointwise or additively, then 
the kernel of K~^ has the same type of decay, respectively. 



1.2 Improving the standard OFDM communications scheme 

In the second part of this dissertation, beginning with Chapter 5, we study orthogo- 
nal frequency-division multiplexing and introduce some extensions of the traditional 
OFDM scheme. In this paper, we develop a comprehensive framework for the design 
of OFDM systems by introducing new pulse-shaping and lattice design methods. 

Transmission over wireless channels can lead to two kinds of dispersion, namely inter- 
symbol interference (ISI) and interchannel interference (ICI). ISI is caused (in part) 
by time dispersion due to multipath propagation and ICI results from frequency dis- 
persion due to the Doppler effect. Two considerations affect the dispersion robustness 
of a given OFDM system, namely the time-frequency (TF) localization of the pulse 
shape, and the time-frequency lattice from which the data are transmitted. 

In Section 6.4, we discuss pulse shaping. We begin with the Gaussian due to its 
optimal time-frequency localization, and orthogonalize via a clever method first dis- 
covered by the chemist Lowdin. This yields a pulse which retains as much of the 
excellent TF localization of the Gaussian as possible under orthogonalization, a fact 
which we prove by appealing to Gabor analysis techniques. Numerical experiments, 
the results of which are presented in section 6.11, illustrate the advantage of the 
Lowdin-orthogonalized Gaussian over the standard square-wave in the presence of 
the Doppler effect. 

OFDM is a special case of Biorthogonal Frequency-Division Multiplexing or BFDM. 
In BFDM in general, one uses a different set of pulses for demodulation than the 
pulses used for transmission. In Section 6.9, we briefly discuss BFDM, in which 
one sacrifices the orthogonality requirement but in return can be afforded improved 
time-frequency localization of the transmitter pulse (one can use a Gaussian, in fact). 
This may improve robustness to multipath and Doppler dispersion, since the disper- 
sion should affect only the transmitter pulse. We then present an extension of the 
BFDM framework which we call extended BFDM, or EBFDM, which is a continuous 
interpolation between OFDM and BFDM. Of all BFDM systems, OFDM, which uses 
the same pulses for transmission and demodulation, is optimally robust against addi- 
tive white gaussian noise (AWGN). EBFDM should have the attribute of permitting 
the communications system designer to fine-tune the balance between robustness to 
ICI/ISI and robustness to AWGN. Inconclusive numerical studies are presented in 
Section 6.12. 

In Chapter 7, we introduce Lattice-OFDM (LOFDM) in which the lattice parame- 
ters are based on sphere-packing considerations (thus maximizing spectral efficiency, 
or alternatively, maximizing robustness to ICI/ISI) and optimally adapted to chan- 
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nel conditions described by he channel scattering function. We complete the pulse- 
shaping by applying a sequence of unitary operators, derived from the metaplectic 
representation of the symplectic group (see Appendix B.3), which optimally adapt the 
pulse to the lattice and channel. Our framework affords the OFDM system designer 
an optimal approach to combatting dispersion. Numerical results are presented to 
demonstrate the superiority of LOFDM in doubly dispersive channels. 



Chapter 2 

Preliminaries and Notation 



2.1 Basics 

Let Z, Q, M, and C denote the integers, the rational numbers, the real numbers and 
the complex numbers, respectively, and set i = \^—l. Let y = a + bi E C; we denote 
by y the complex conjugate y = a — hi. 

Denote by Vt the set of all functions defined on M?'^. 

"There exists" will be deonted by the symbol 3, and "for all" by V. The set with no 
elements, called the empty set, is denoted by 0. 

A vector space or linear space X over a scalar field M is a set of elements on which 
are defined operations of elementwise addition (+) and (commutative) scalar multi- 
plication with the following properties: 

(a) the operation + is commutative, i.e., for all x,y E X, we have x + y = y + x; 

(b) the operation + is associative, i.e., for all x,y,z G X we have {x + y) + z = 
x + {y + z); 

(c) there is a zero vector such that for all a: G X we have x + = x; 

(d) for all X G there exists a unique element —x such that x + {—x) = 0; 

(e) since M is a field, it contains the element 1, for which Ix = x for all a; G X; 

(f) for all x,y E X and A,/i G M we have (A + fi)x = Xx + fix, {Xfi)x = X{fix), and 
X{x + y) = Xx + Xy. 

A norm on a linear space X is a function || ■ || ^ M satisfying the properties 

(a) ||a;|| > for all x G X; 

(b) ||Aa;|| = |A|||a;|| for all A G C and a; G X; 
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(c) ||a; + y\\ < ||a;|| + ||i/|| for all x,y & X; 

(d) ||a;|| = <^ a; = 0. 

If a function / is absolutely p*^-power integrable, we denote the p-norm of / by 

where 

1 

p 



VI 



\f\x)\^dxj . (2.1) 

Note that for two functions /i and /2 we can have ||/i — /2||p = even if there exists x 
such that fi{x) ^ f2{x). Thus is generated an equivalence class structure. The space 
of equivalence classes of functions / for which ||/||p < c)o is denoted L (M'^).^ 

The essential supremum of a nonnegative function /, denoted ess sup f{x), is defined 

by 

ess sup f{x) = inf {M I m{x \ f\x) > M} = O} (2.2) 

where m{T) denotes the measure of a set F. Again, the essential supremum does 
not capture pointwise differences between two functions, and generates equivalence 
classes. The space of equivalence classes of functions / for which 

||/||oo:=ess sup|/(a;)| < oo (2.3) 

is denoted L {W^). Not all interesting functions are in L {W^) or L (M'^). The 
space of infinitely differentiable functions on M^ is denoted C°^(]R'^), and Cq := 
{/ I / continuous and /(x) — ;► as ||a;|| -^ oo}. 

2.2 Banach and Hilbert spaces 

A normed linear space (X, || ■ ||) is a linear space X equipped with a norm || ■ ||. A 
Cauchy sequence {xn} is a sequence in which for every e > there exists an A^ G N 
such that for all n,m > N, we have ||a:„ — Xm\\ < e. A linear space X is complete 
if every Cauchy sequence in X converges to an element in X. A Banach space is a 
complete normed linear space. Familiar examples of Banach spaces are the L (M^), 
with 1 < p < oo. 

Lemma 2.2.1 Let Bi and B2 are Banach spaces with the same metric d. Then Bir\B2 
is a Banach space. 

Proof. Let i3i, B2 be Banach spaces with metric d such that i3i fl ^^2 7^ 0- We need 
to show that every Cauchy sequence in Bi fl B2 converges to an element a; in i3i fl i32. 
Let {xnjneN C Bi n B2- Then since {x„} C Bi, {xn} converges to an element x in 

-"^We shall follow the usual convention in denoting equivalence classes by /, but referring to / as a "function." 



Bi. By virtue of its membership in B2, {xn} converges to an element x' in 82- If 
d{x,x') = e > 0, then 3e' and A^ such that either d{xp^,x) > e' or d{xiy,x) > e' for 
all n > N, contradicting the Cauchy property of {a;„}. Thus x' = x. We conclude 
that {xn} converges to a unique element in Bi ni32, and so Bi nB2 is complete in the 
metric d, and as such is a Banach space with norm inherited from Bi and B2- Q 

An inner product on a linear space X is a map (■, ■) : X x X ^ C with the properties 
that for all x,y,z ^ X and A, /i G C, we have 

(a) {x , \y + jiz) = \{x , y) + ^{x , z) 

(b) (y, x) = (x, y) 

(c) (x, x) >0 

(d) (x, x) = <(^ X = 0. 

We call a linear space equipped with an inner product an inner product space. If X 
is an inner product space, we can define a norm on X by ||x|| = a/(x, x). A Hilbert 
space is a inner product space which is complete under the norm derived from the 
inner product. Famous examples of Hilbert spaces are L^(]R'^) and £^(Z'^), the space 
of square-summable sequences of complex rf-tuples. 

2.3 Algebras 

An algebra is a vector space A on which the following properties are satisfied: 

a{bc) = {ab)c (2.4) 

a{b + c) = ab + ac (2.5) 

{a + b)c = ac + bc (2.6) 

X{ab) = {Xa)b = a{Xb) (2.7) 

for each a,b,c & A and scalar A. 

An involutive algebra, is an algebra A on which there exists a unary operation (*) 
such that 

(2.8) 

(2.9) 

(2.10) 

(2.11) 

for all a E A. 



a 


= a 


a + by 


= a* + b* 


{Xay 


= X*a* 


(ab)* 


= b*a*, 
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We call A a unital algebra if A contains the multiplicative identity. 

A normed algebra is an algebra A on which is defined a norm || ■ || which is submulti- 
plicative, i.e., 

\\ab\\ < \\a\\\\b\\ 

for all a,b E A. 

A Banach algebra is a normed algebra which is also a Banach space, i.e., complete in 
the norm 11 ■ 11. 



2.4 Fourier considerations 

The Fourier transform of f E L^ {M.'^) denoted by / is 

/(.;) = / /(t)e-2™-*rft. (2.12) 

We also write J-'f instead of /. If / G L^(]R'^), then we have the inversion formula 

fit) = [ /(c.)e-2™-Vc.. (2.13) 



The Fourier transform can be extended to a unitary operator on L^(]R^); from here 
on the notation JF shall represent this unitary operator where applicable. 

Fix a function g E L (M'^), q> 1, (^ 7^ 0, and let p = -^. Then the short-time Fourier 
transform (STFT) of a function f E L (M.'^) is given by 



Vgf{x,uj)= f{t)g{t-x)e'^^'''-'dt, x^ueR'^. (2.14) 

The ambiguity function of / G L'^{W^) is defined to be 

Af{x,uj) = f f{t + ^)/(t-^)e-^-^-*rft, (2.15) 



where f{t) denotes the conjugate of f{t). Let p > I , q = -^. The cross-ambiguity 
function of / G L\R'^),g G L^(M'^) is 



Aif,g)ix,co)= I /(t + |)^(t-|)e-2--*rft = e-^--^\/,/(a:,c.). (2.16) 



10 

2.5 Tempered distributions 

Let Z_|_ = {n G Z I n > 0} denote the nonnegative integers. A multi-index a = 
(cti, . . . , ttrf) G Z^ is a (i-tuple of nonnegative integers. Let x = (xi, . . . , Xd) G M'^ and 
a G Zi, and define 



D' 



"^1 / ;5 \ "^d 



9 \"^ / 9 



dxi J V Sxji 



d 



a 



n 



N = V^? + --- + ^^ (2.17) 

We will use the notation in (2.17) for the finite-dimensional 2-norm throughout the 
sequel. The Schwartz space S consists of all ip G C°^(]R^) such that sup |a;°D'^(/9(a;)| 



is finite for every pair of multi-indices a,/? G Z^. Elements of S thus decay faster 
than the inverse of any polynomial as |a;| -^ oo. 

The space of continuous linear functionals on S, or the topological dual space of S, is 
denoted S*. Elements of S* are called tempered distributions. The flagship example 
of a tempered distribution is the Dirac delta "function" defined by 

5i^) = vp(O). 
The most intuitive way to view 5 is as a measure on S, whereby 

ip{x)6{x) dx = (p{0). 



We can extend the domain of the point measure 6{x) to any function defined on M'^ 
(including the constant function 1) enabling us to say f^^^ 6{x) dx = 1. Note that in the 
sense of Lebesgue, f^^ 6{x) dx = 0. The Fourier transform of a tempered distribution 
ip can be defined by (p{f) = V'(/) where / = J^~^{f), f E S. Thus 

S{f) = f S{x)f{x) dx = /(O) = / -^/(cc^) doo (2.18) 

JRd jRd v27r 

which implies that, in the sense of measure, 6{uj) can be thought of as the constant 
^4=. Thus 6{u) is a multiple of the usual Lebesgue measure, and S{x) can be thought 



of as the inverse Fourier transform of the Lebesgue measure, up to a constant. 
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2.6 Integral operators 

Let /i be a member of some function space X on M"'. An integral operator on X is an 
operator K defined by 

K{h){x):= [ k{x,y)h{y)dy (2.19) 



when Kx G X, i.e. K : X -^ X. We call k{x,y), x, y G M"* the kernel of -ftT. In this 
paper we shall remove the distinction between K and k since k defines K uniquely, 
and will refer to the kernel and the operator by the same symbol. That is, in the 
sequel k{h) = K{h). 

Let / : M'^ ^ M be any function. The Laurent function on M?'^ induced by / is 
fhix^y) = f{x — y). We say that / generates fi. This is the bi-infinite analog of 
Toeplitz functions on M.'^, and the idea of a Laurent function is of course well-defined 
in discrete settings as well. The Laurent function fi can be viewed as the kernel of 
a convolution operator, provided that the convolution is defined. Let g he & function 
defined on M'^. Then the convolution of / with g is defined to be 

f*9it)-= [ f{x)g{t-x)dx= [ f{t-x)g{x)dx, (2.20) 

provided the integral in (2.20) is defined. If we view f^ as an operator, we have 

f*g = fdg)- 

2.7 Continuous analog of matrix multiplication 

We write explicitly the integral operator formed by composition of integral operators. 
Let I < p,q < oo, - + - = 1, f = f(t,y), g = g{x,t), and define the restriction operator 
R as R^2d{h){x,y) := h\^2d{x,y). We further demand that 

(1) each "row" in the matrix A := R^2d{f){t,y) is absolutely p-summable, and 

(2) each "column" in the matrix B := R^2d{g)[x,t) is absolutely g-summable 

(making Rj2d{f) and Rj2d{g) bounded linear operators on (i {!/) and (i (Z'^), respec- 
tively). For A, let the row index (y-axis for /) be j and the column index (t-axis for 
/) be £, and for B let the row index (t-axis for g) be £ and the column index (x-axis 
for g) be k. We have 

{AB)^k = J2A,,B,k. (2.21) 

It is clear that the analogous continuous multiplication, defined for the moment on 
the space of all functions g : M^'^ -^ M (and for the moment ignoring convergence 
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considerations), denoted as (^(M^*^) which respects axes and orientation must be 

{f^g){x,y):= f f{t,y)g{x,t) dt, (2.22) 



where we shall assume that x = Qii g E VL{W^). To visualize the difference in indexing 
between this operation and matrix multiplication, consider g G r2(]R). In this case, 
the convention should be that g is visualized vertically, with positive y-values at the 
top and negative y-values at the bottom. Of course the operation -k satisfies 

U^9)^h = f^{g^h), (2.23) 

when both f -k g and g -kh exist. A proof of associativity is given in Appendix A.l. 

We shall later show that certain classes of tempered distributions in S* act on the 
Hilbert space L'^iW'-). In this case, we have, for hi, h2 E L^(]R^), 



(/ii,/M = / / h{y)f{t,y)h2{t)dtdy 



hiy)fiit,y)dyh2it)dt 



fiy,t)hi{y)dyh2it)dt 

= (r/^i, M 

which implies that /* = (f)^- Thus, the adjoint operation in any subspace of ^2 which 
acts on L^(M'^) is given by refiection about the line t = y followed by conjugation, 

i-e., 

nt,y) = f{y,t). (2.24) 

In addition, we have (abusing notation slightly) (/ -k g)* = g* -k /*, from 

if^gT = ( [ fit,y)gix,t)dt^ 

T 



fit,y)g{x,t)dt ^ 
fit,y)g{x,t)dt 



git,x)fiy,t)dt. (2.25) 



2.8 •-invertibility 



The left unit for the multiplication operation -k is 6{t — y); the right unit is S{x — t); 
clearly these are equal. But -k need not be a group operation. By illustrating the action 
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of a convolution operator /^ on a function in L^(]R'^) which "oscillates to death" in 
the sense of Lieb and Loss [37], we will show that f^ is not boundedly invertible. 
Consider the Gaussian /(t) = e^'^'*' , t G M*^, and form the Laurent operator fi on 
L2(M°') with 

fl{t,y)=e--\y-'r (2.26) 

Note that /l is not Hilbert-Schmidt (the kernel fi ^ L^(]R^'^); Hilbert-Schmidt opera- 
tors are compact and hence not invertible). However, we do have inf ||/l^5' II2 = 0, 

which precludes invertibility by a bounded linear operator. Indeed, consider g^^ = 

{2i)f{t)e^^'^-* which has ||^^||2 = l,Vcj G W^. We have 

= J^(e-"l^-*l'e-"l*l')(i/,cc;) 

which is just the STFT of / with / as its own window. Since ||V//||j^2(-]R2d) = 
Ii2(]jd) < 00, we infer that the slices \\Vff{y,Lj = a;o)||L2(iRd) ^ as \uo\ -^ 00 



which implies W/l'*^ 9u)\\L'2{R'i) ^ as |cj| ^ 00 (this result is related to the Riemann- 
Lebesgue lemma). Proposition 4.8 in [16] states that if T is a bounded operator on 
a Hilbert space, then T is (boundedly) invertible on the Hilbert space 2y^(]R) if and 
only if T is bounded below and has dense range; thus we have shown that fi is not 
invertible. 

In the sequel we will drop the symbol -k and simply denote f -^ g ^y fg when the 
context of the multiplication is clear. 

We next consider the existence of invertible integral operators, specifically invertible 
convolution operators. Let /, g be functions in function spaces yet to be determined. 
We define the convolution of / with g to be 

{f*g){t)= I f{x)g{t-x)dx. (2.27) 

We wish to invert this operation on g. Young's inequality says that ii f,g G L^{M.'^) 
then f * g & L^{W^). In this case we have 

T^g = f-g (2.28) 

and thus ^ 

g = ^-^ (^) = :f-' (t^) = :F-\g) (2.29) 

showing that the convolution operation on g is invertible (since / G Co(]R'^), we can 
use a limiting process at the points where / = 0). 
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But what about the operator /~^ in general, operating on functions not necessar- 
ily formed by convolution? When is it a bounded operator on, say, L^(]R^)?^ By 
inspection of (2.29) we have 

Ifinf |/H|>0, 

iM < J(^ = Chiu) (2.30) 

ficu) inf/(c.) 

which is in the domain of JF~^ if /;, G L^(]R^). Obviously if lim f{uj) = or 
lim f(uj) = 0, then 4^1 is not bounded, hence not integrable and not in the domain 

of T^^ . It is also possible that, for a given /i, -^^ E L^(R.^) with f{uj) = at a count- 

able number of points Un- Consider, e.g., an / which is different iable, and whose slope 
is sufficiently large, and sufficiently increasing as n ^ ±oo, at the zero crossings of 
f{uj) so as to ensure that -^^ G L^(R.'^). For example, suppose that (loosely speaking) 

away from the zero crossings, / = 1, and near the zero crossings j^ had the same 

type of vertical decay behavior as {x — Un)^^^ (for any e > 0, /'^"_ ''{x — Uny^^ dx 
exists and as n ^ oo is decreasing sufficiently rapidly so that the sum of all the 
integrals about the countable number of singularities is finite). But no such / will 

satisfy i^ G L^R"^) for all h G L^R"^). 

So what kind of function is / with inf f{uj) > 0? The Riemann-Lebesgue lemma states 

that h G Li(]R'^) ^he Co(]R'^). Thus neither / nor /can be an /.^(M"^) function. 
Could / be a member of another L (M^) space? Certainly not. Let C = inf f{uj). 

LU 

Then J-'~^f = -S= S + /2, where /2 may be in some L (]R°') space but 5 ^ L (M"') for 
any p since 6 = 6 is not a function in the Lebesgue sense (see Section 2.5). 

Thus for a convolution operator to be invertible, it is necessary that the kernel of the 
operator possess as a summand a Dirac delta distribution. I.e., 

fL{x,y) = f{x-y) with f = f + S. (2.31) 

We will call / the generator oi the operator f^. In the discrete setting, the analogous 
condition is not sufficient; this is called the Wiener-Pitt phenomenon. For any locally 

^We could ask when /^^ is a bounded operator on the Hilbert space L^(R'') but for the purposes of the present 
argument, we will not complicate matters by appealing to the limit arguments necessary for the treatment of J^ and 
jc-l in L2 
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compact abelian group G, define the Wiener algebra W{G) as 

W{G) = {geL\R^)\geL\R^)}. (2.32) 

In our setting, it is known that if the nondistributional part / of the generator / of the 
Laurent operator is in the Wiener algebra W{G), then the Wiener-Pitt phenomenon 
does not occur and condition (2.31) is sufficient as well as necessary. We assume in 
the sequel that if an operator g^ under consideration is a convolution operator, its 
generator g has nondistributional part g G W^(]R'^) C L-'^(]R'^). We can then safely say 
that gi acts on /.^(M^). 

2.9 Spaces of distributions 

Recall that for a convolution operator g^ to be invertible, we must have gL{x,y) = 
g{x — y) + 6{x — y), where (7 : M*^ ^ M. We extend that idea now to the kernel of any 
integral operator. Formally adjoin to Q the unit element e := 6{x — y) to yield the 
class Q, whose elements are of the form f = f + /ue, f E Q , fi eR. We shall consider 
the class Qi, whose elements are of the form / = f + fie + X^feli CkS{xk,yk), f ^ 
n , /i , Cfc G M, and where each hyperplane defined by constant x or y supports only a 
finite number of deltas. 
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Chapter 3 

Classes of Integral Operators 
Classified by Off-Diagonal Decay of 
the Kernel 



In this chapter we will show that if the kernel of an integral operator possesses 
subexponential decay, whether measured by weighted "pointwise" essential supre- 
mum boundedness, or by essential supremum boundedness of the weighted partial 
integral, then the kernel of the inverse integral operator possesses the same type of 
decay, measured in the same manner as the original kernel. 

3.1 Weight functions 

We define the relevant classes of distributions under consideration by their decay 
relative to weight functions. We measure the decay of an integral operator kernel in 
two ways: by boundedness of its weighted integral or pointwise essential supremum 
boundedness against weights. Weights are used in order to quantify the convergence 
properties. In addition, the quantitative nature of measurement against particular 
weights will be of some advantage when deriving numerical methods and assessing 
their performance. 

Definition 3.1.1 An admissible weight v is (a) subexponential: 

v{x) = e''''^'' with K concave and k{0) = 0, and (3.1) 

(b) satisfies the Gel'fand-Raikov-Shilov (GRS)-condition (cf. [20]) 

lim t;(nx)^/" = 1 Vx G R'^. (3.2) 

n— >oo 

We shall discuss only admissible weights, which lie in L^^^{M.'^) , 1 < p < oo, in what 
follows. An example is Ts{x) = (1 + |a;|)*, s > 2d. But to further clarify what sort 
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of function is an admissible weight, we state some of the imphed properties of such 
weights. Note that since k is concave and /t(0) = 0, f is submultiplicative, i.e., 

v{x + y) < v{x)v{y). (3.3) 

Although K{t),t G M is concave v{x) = e'*'-''^'-' is neither log-concave nor log-convex. 
In fact, if X G M^, then every point on the graph of v{x) resembles a saddle point (see 
Figure A.l in Appendix A. 3). As a consequence of the concavity of v and the fact 
that K, is nonnegative, k must be nondecreasing (increasing if v is strictly concave). 
Finally, we note that since k, is concave, v must be continuous. For a more detailed 
analysis of each of these properties and an illustration of the GRS-condition, see Ap- 
pendices A. 3 and A. 4, respectively. 



3.2 The distribution classes jCI and Cy 

We define the relevant classes of distributions. 

Definition 3.2.1 The class Cl consists of all elements f E Qi such that 

Sx,vif) '■= ess sup / \f{x,y)\v{x-y) dy < oo, and 

Sy,vif) '■= ess sup / \f{x,y)\v{x — y)dx<oo. 

Define a norm on Cl as 

WfWd := max{5.,,(/) , SyM)} + |/i| = ||/|Ui + |/i|. (3.4) 

Write £\ Sx{f), Sy{f) in the case of f = 1. Note that / G £^ is quite similar to a Hille- 
Tamarkin operator (cf. e.g., [35]), the difference being that integrability of members 
of Cl is measured against weights, and that elements of Cl may contain the identity 
and/or point measures. Alternatively, Cl can be (nearly) viewed in the context of 
the L-^(]R^^) spaces of mixed-norm (cf.[4]).^ Obviously the above relationships are 
stronger when v = 1. 

Definition 3.2.2 The class C^ consists of all elements f E Qi such that 

ess sup\f{x,y)\v{x—y)<oo. 



Define a norm on C^ as 

\c, :=ess sup\f{x,y)\v{x-y) + |/i| = ||/||£„ + |/i|. (3.5) 
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Figure 3.1: Depictions of truncated elements of £^ or £„, derived from compact and noncompact 
operators 



In terms of Hadamard multiplication, if we let vl{x, y) = v{x — y), we have 
£,„=</ G f2i ess sup / (/ vl) dx<oo and ess sup {f Q vl) dy< oo 



and 



C 



fen, 



ess sup (/ Vl) (x, y) < oo ^ . 



The inclusion of a infinite number of delta-distributions, with only a finite number 
supported on any given hyperplane defined by constant x or y, in members of Cl and 
£y is permitted since the deltas make only finite contributions to the integrals on 
sets of measure zero, thus preserving the equivalence class structure of L^(]R'^) and 
inducing equivalence class structures on both Cl and £„. We shall adopt the same 
convention as one does for the L (M'^)-spaces, namely that although / in, say, Cl 
represents an equivalence class, we will take liberties with the terminology and say 
"/ is a operator in £j." 



It is also interesting to note that neither Cl nor Cy are separable (Appendix A. 6). 



3.2.1 The unital Banach algebra Cl 

Lemma 3.2.3 // v is admissible, then Cl is a unital Banach *-algebra of bounded 
operators acting on //^(M*^). 



^Note the capital P, generally signifying P = {pi},i = 1, . . . n with '}27=iP' ~ ^' ™ °^^ case n = 2. 
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Proof. Considering first the unweighted case, if / G £^, then / is simultane- 
ously bounded on L {W^) {f^^ \f{x,y)\ dx is bounded) and L^(]R'^) {f^^ \f{x,y)\ dy 
is bounded). Thus / is bounded on L (M*^), 1 < p < cxo, and in particular L'^(W^), 
by Schur's test (see Appendix A. 2). So C^ is a set of operators bounded on L'^iM.'^), 
but is it an algebra? Yes, the algebra axioms ( (2.4) - (2.7)) are easily verified. In 
fact C^ is a *-algebra. Also, C^ is a Banach space, since C^ coincides with L{^ of 
[22], which was shown to be a Banach space in, e.g., [4]. It remains to show that 
C^ is a normed algebra; i.e., that || ■ H^i is submultiphcative. Let f,g E £\ with 
f = f + /iie, g = g + /i2e. The definition of || ■ H^i yields 



WfgWc^ = max|||/^||op,oo, ||/^||op,i \ + |^i||/i2| 
and we have 

||/^||op,oo = sup \\fgh\\oo 

l|ft||oo = l 

\\m\ 

= sup 



II lloo 

llMlloo \\~9hl 

SUp -T-TT, Tr^TTT 

heL'^im i|/i||oo wg'T'Wc 



Wf-ghWoo \\~gh\\c 

oLlL) 

heL° 



^°°(Rd) ll^^lloo \\h\\c 



< 



WfghWoo \\gh\\c 

sup — sup 



ll^/^ll 
11 

h 



oo heL'^ 



II Alloc ,1-11 
< sup — — I ||5'||op,oo 



oo 



|op,oo II y II op, oo 

Thus the cx)-norm is submultiphcative. The 1-norm is also submultiphcative, proven 
in the same manner. And because 



max 



op,i 



SOo||fi'||op,oo , ||/||op,l||fi'||op,l ( + |/^l||/^2 



< max<^ ||j||op,oo, ||j||op,l } + IfJ'l 



■ max<^ ||5'||op,oo , ||5'||op,i > + |/i2| , 

we conclude that || ■ H^i is submultiphcative, which establishes that C^ is a normed 
algebra, and thus a Banach algebra (with unit by contruction). 
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For the weighted case, we begin by showing that it is reasonable for / G £^ to act on 
L^(M'^), by first demonstrating that / acts on L'^{M.'^), then extending the conclusion 
to all of Cl- We can no longer use Schur's Test. Let h G L^(]R'^), with \\h\\p = 1 (it 
is easy to see why no generality is lost under this assumption). Then the action of 
f G Cl on h E L {M.'^), 1 < p < oo, is given by 



(my) 



~f{t,y)h{t)dt. 



(3.6) 



Let Qhit^y) be the Laurent operator over M^'^ formed from translations of 
9a{y) = f { arg<^maxiess sup / \f{t,x)\\h{t)\dt } \ , y) . 

I.e., we form the Laurent operator gL{t,y) := ga{y — t) from the slice of / which 

acts on h to yield the largest ess sup Jj^^j \f(t,y)\\h(t)\ dt. It is clear that ql inherits 

yes.'' 

membership in Cl from /, and we can conclude that Qa G -/^^(M"*). Then we have, for 
1 < p < oo. 



IIAII^ 



< 



\miyWdy 



f{t,y)h{t)dt 



dy 



\f{t,y)\\h{t)\dt] dy 



< 



\gL{t,y)\\h{t)\dt] dy 



(3.7) 



\gaiy-t)\\hit)\dt] dy 



\g,\ * \h\ (y) r dy 



= m9a\* 

< \\9a\mh\\; 



P \\h\\P < oo 



(3.8) 
(3.9) 



where (3.7) follows from the definition of g, (3.8) follows from Young's inequality 
with r = p and q = 1, and (3.9) follows from the definition of ga{y) and the fact that 
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v{k) > 1. Extending the analysis to general f = f + ^e E C\ where /i 7^ we get 

wmi = \\if+Hh\\i 

< \\fh\\l + \f^nh\\l 

< ll/T^jll/^ll^ + 1/^111/^11^ 

< WfWimi (recall p>l). (3.10) 

Thus for all h G L^{R^), we have \\fh\\p < M\\h\\p and we conclude that any f G Cl 
is a bounded linear operator on L (M.'^) (and in particular L^(]R'^)). 

For Cl, the algebra axioms are all easily verified. To show submultiplicativity, we 
again begin by considering elements of £^ which do not include the unit as a summand. 
We have 



\fg\{x,y)v{x-y)dx 



< 



< 



f{t,y)g{x,t)dt 



v{x — y) dx 



\f{t,y)\ \g{x,t)\v{t-y)v{x-t)dtdx 
\fit,y)\ \g{x,t)\v{t - y) v{x-t) dx dt 



< / \f{t,y)\v{t-y) ( \~g{x,t)\v{x-t)dx] dt 
< / less sup / \g{x,t)\v{x — t) dxj \f{t,y)\v{t—y)dt 

< 



£1 9 £1 

*-'ii II .^ II *-'i: 



and likewise for 



\fg\{x,y)v{x-y)dy. 



Taking appropriate suprema yields H/f?!!/:! < ||/||/:i lls'lUi- Extending the analysis to 
the treatment of elements f = f + /iic , g = g + /i2e ^ ^v, we have 

WfgWci = ||(/ + /iie)(^ + ^2e)|Ui 

= ||/^ + Ati^e + /i2e/ + /ii/i2e||£i 

< WfgWcl + l/^llll^lki + I/^2|||/||£J + |/il||/i2| 

< ll/IUill^Iki + l/^illl^lki + l/^2|||/||/;j + |/ii||/i2| 

= i\\f\\cl + \f^l\)m\cl + \f^2\) 

IclhWcl (3.11) 



showing that || ■ H^i submultiplicative, and we can conclude that Cl is a Banach 
algebra. □ 
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3.2.2 The unital algebra £„ 

Definition 3.2.4 A function g is called subconvolutive if g satisfies g^^*g^^ < Cg^^, 
for some C > 0. See Appendix A. 5 for an elementary example of an admissible weight 
which is subconvolutive. 

Lemma 3.2.5 If J^aV^^{x) dx < oo, then C^ C C^. If in additionv is subconvolutive, 
then Cu is a unital *- algebra of bounded operators acting on I?'(W^). 
Proof. As before, let / = / + /ie, f & C^, with v assumed subconvolutive. We first 
show Cy C C^. Since 

\f{x,y)\v{x-y) < ess sup \f{x',y')\v{x'-y') 



almost everywhere, we conclude that 

\f{x,y)\<\\f\\c.v~\x-y), (3.12) 

except possibly on a set of measure zero. Then if J^v~^{x) dx < oo, we must have 



ess sup 



\f{x,y)\dx 



< 



ess sup 



*—v 



*—v 



\f{x,y) + i2e{x,y)\dx 

v~^(x — y) dx + |/i| 
v~^{x) dx + |/i| < oo. 



(3.13) 
(3.14) 



and similarly for ess sup {/j^ \f{x, y)\ dy] . Thus Cv C C^. Next, we show that C^ is an 

xeR'i 

algebra of bounded operators on L (R'^). We begin by demonstrating that arbitrary 
/ G £„ in fact acts on L'^{W^). Let h G L'^(M.'^). Recall the (formal) action of / G £t, 
on h is given by 



ifhM 



Then 



KfhM^dy 



f{t,y)h{t)dt. 



f{t,y)h{t)dt dy 



< 



< 



p 
p 



v-\t - y)h{t) dt 



dy 



\v-'*h\\l 



,-i||P 



\h\\P < oo 



(3.15) 

(3.16) 
(3.17) 



where (3.15) follows from the estimate (3.12), (3.16) is a consequence of the symmetry 
of V, and (3.17) follows from Young's inequality with r = q and p = 1. Thus for all 
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h E L (W^), \\fh\\p < oo and we conclude that any f E C^ induces a bounded linear 



operator on L 



(in particular for p = 2). To show that Cy is an algebra, we need 



to show that the class £„ is closed under multiplication. Without loss of generality 
we can assume the submultiplicativity constant C > 1. We have 



\{m^,y)\ 



< 



f~'V W tJ W f~'V 



f{t,y)g{x,t)dt 



v^\t-y)v-\x-t)dt 



< \\f\\cM\c.Cv-\x-y) 

^\{r9){x,yMx-y) < cWfWcSgWc 
^ WMc. < C\\f\\cXg\\c., 



(3.18) 
(3.19) 

(3.20) 



where (3.19) follows from the subconvolutivity off, and (3.18) and (3.20) hold except 
possibly on a set of measure zero. Invoking (3.20) in a calculation nearly identical to 
(3.11) yields 

WfgWc. 



< \\mci + w\mci + \^2\\\f\\c. + \f^i\\f^2\ 

< C ||/||/;i||^||/;i + l/iilll^lUj + |/i2|||/||/:„ + |/il||^2| 

< C ll/IUill^lUi + C \fii\\\g\\ci + C l/^2|||/IU„ + C I/X1II/X2I 

= C(||/|U„ + |/ii|)(||^|U, + |/X2|) 

= CWfWcMlc. (3.21) 



Thus fg satisfies the definition of C^. The class C^ possesses a unit by contruction, 
and thus Cy is a unital algebra of bounded linear operators on L^(]R°'). □ 

We see that under the norm || ■ H^^ and the operator multiplication fg, C^ is not 
strictly submultiplicative, so we seek a different norm under which C^ might satisfy 
submultiplicativity. A reasonable candidate might be || ■ H'^^ defined on C^ by 



*—v 



ess sup \\fh\\c. 



(3.22) 



hGCv , II h\ 



Cv~- 



which is the continuous analog to the norm || ■ ||^^ mentioned briefly in [23]. Unfortu- 
nately, the two norms || • \\c^ and || • \\'^^ are not equivalent. In [23] the corresponding 
norms || ■ ||^ and || ■ ||_4^ are equivalent but the relationship breaks down in the con- 



tinuous case; consider the sequence {/„}neN where /„ 
||/n||£„ = lVn, but lim||/„||'^^^0. 



[X 






'x). We have 



24 



Transforming L^ into a Banach space 

We define a new norm and space, and show that, under broad circumstances, under 
the new norm, the new space is equal to £„. 

Fix s > 0, and let u be an admissible weight on M'^. Let T^ix) = (1 + |a;|)*, (also 
admissible) and set 

v{x) = u{x)rs{x). (3.23) 

Note that 

Tsix + y) < K{Tsix) + Tsiy)) (3.24) 

for some k. To prove this we consider two cases (in each we assume x,y > 0): 



X > y 



X < y 



1 + , + ,). ^(1 + 2^^ ^^^ 



{i + xy + {i + yy {i + xY 

{l + x + yY ^{l + 2y) 



< 



[i + xy + {i + yy - (l + yY 

note that k, > 1. We also have 

v{x+y) = u{x+y)Ts{x+y) 

< Ku{x)u{y){Ts{x) + Ts{y)) 
= K[v{x)u{y) + u{x)v{y)]. 

Define the Banach space B^^s '■= ^l. f^ ^v with norm 

U. + l/^l- 



<2' 



Bu,s = HJWCI 



(3.25) 



(3.26) 



Lemma 3.2.6 (a) The class B^^s with norm (3.26) is an involutive unital Banach 

algebra. 

(h) If s > d, then Cy C £^ and thus Bu^s = '^v with equivalent norms. 

Proof, (a) Suppose f,g E Bu,s- We have 



wmi 



ess sup 



fit,y)gix,t)dt v{x-y) 

< ess sup / K\f{t,y)\v{t — y)\g{x.,t)\u{x — t) dt 
+ ess sup / K\f{t,y)\u{t — y)\g{x,t)\v{x — t)dt 

x,y<^R'i Jr'' 

< «;||/IU„ess sup / \g{x,t)\u{x-t)dt 
+ '«ll^lk.ess sup / \f{t,y)\u{t-y)dt 
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< '^(ll/lkJI^Iki + ll/lkll^lUJ- (3-27) 

Thus it follows that 

< i^{\\f\\cj9\\ci + WfWciMcJ 

+ |/i2|||/||£, + |/^l|||^|k, + |/il||/i2| 

< <\\f\\cM\ci + \\f\\ci\\9\U- (3.28) 

(We actually do not need (3.28) to show that || ■ ||g„ ,, is submultiplicative, but we 
shall need it in Section 3.3.4). We also find 

\\{f9)\\c. < ||/^|k„ + |/i2|||/|U„ + |;Ui|||^|U, 

< <\\ncA~9\\ci + \\f\\ci\\9\U 

+ |/i2|||/||i:. + |/ii|||^IU„. (3.29) 

To establish that || ■ \\bus i^ submultiplicative, we require that 

||(7^)IUi<||/IUill^lki+/^2||/lki+/ii||^IUi, (3.30) 

which is easily proven using the definitions. Using (3.29) and (3.30), we have 

\\f9\W,. = 4(f9)\\ci + \\G9)\\c. + \f^i\\f^2\ 

< '^dl/llall^lki +^2||/||£i +^l||^|Ui) 



t^[\\J\\cJ\9\\ci + ll/lkill^lkJ+^2||J||/:, +/ii||^||/:, 

|/ii||/i2|. (3.31) 



BuJ\9\\bu,s = '^'^WfWcihWcl 



But 



+/^(ll/IUill^lk„ + l^2|||/IUi) 

+ ||/||/:J|^IU„ + |/i2|||/IU„ 

+ lAti|||^||/:. + |/Wi||At2|. (3.32) 

And since 11 • 11 ri and 11 ■ llr are norms and thus nonnegative, and /t > 1, we imme- 
diately obtain ||/^||b„,, < H/b^.J^He^ .. 
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The adjoint operation is an involution on Bu,s- The intersection of two Banach spaces 
is again a Banach space as we saw in Section 2.2; thus B^^s with norm (3.26) is an 
involutive Banach algebra. 

(b) li s > d, then t^^{x) is integrable, so 

WfWci = \\f + f^e\\ci 

= ess sup / \f{t,y)\u{t-y)Ts{t-y)T~^{t-y)dt+\fj,\ 

< WfWc. f T:\t- y)dt+ |/i| 

= A\\f\\c^ + U (3.33) 

Now let A' = max{l, A}. Then we have ||/||£i < A' \\f\\cv ^-iid thus the norms on 
Bus and C^ are equivalent: 



\c. < 4f\\ci + ll/IU„ = \\fh.,s < M' + 1)11/11/:., 
and Cy C C^ ^^ Bu^s = ^v Q 

3.3 Symmetry, and decay of inverses 

3.3.1 Symmetry of algebras 

Before we show the main result on the decay of 7»r-inverse kernels, we return briefly to 
the question of algebra symmetry, since the symmetry results shall be proven along 
the way. Recall that o-_4(/) is the spectrum of / in the algebra A, and a{f) is the 
spectrum of the operator / acting on a Hilbert space 7i. In general, (TaW) ¥" ^(/)j if 
you can read French, see reference (9) in [25]. The respective spectral radii are 

PaU) = max{|A| : A G o-^(/)} 
p(/) = max{|A| : \ e a{f)} . 

If / is self-adjoint, then p(/) = ||/||op (see, e.g., [39]). 

Lemma 3.3.1 If A C B, then we have PbU) < PaH)^ V/ G ^. 

Proof. For A G o'Aif) there may exist an element in B\A which inverts XI — f. □ 

A pair of nested algebras ^ C i3 is called a Wiener pair if / G ^ and /^^ G B implies 
f^^ G ^. ^ is often termed inverse closed in B. 
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A Banach algebra A is called symmetric if crA{f*f) C [0, oo) \/ f E A. Equivalently, A 
is symmetric '^ f = f* implies (Tj\^{f) C M. Recall that if / = /* is a bounded linear 
operator on a Hilbert space, then a{f) C M (see, e.g., [39]), so if we can estabhsh for 
a particular A that V / = /* we have o"_4(/) = ct(/), we will have shown that A is 
symmetric. 

3.3.2 Crucial lemmas 

Lemma 3.3.2 (Hulanicki [25]): If A is an involutive Banach algebra with identity 
contained in BiTi), and if for all f = f* E A we have 

PAif) = P{f) = WfWo, (3.34) 

then aM) = ^(/) for all f = f* e A. 

As a consequence of Hulanicki's lemma, we conclude that ^ is a symmetric Banach 
algebra. In addition, A contains the inverses of each of its self-adjoint elements. To 
see this, suppose W f = f* E A,a{f) = o'_4(/). Then if / is not invertible on A, 
we have G o"^(/) and thus e a{f) so / is not invertible on BiTi). Thus by the 
contrapositive we conclude that if / is invertible on B{7i), then / is invertible on A. 
If we can extend this latter conclusion to all of A, then we will have shown that A 
and B{7i) form a Wiener pair. 

Lemma 3.3.3 (Barnes [2]): Suppose that w{x) = (1 + \x\Y for some 6 G (0, 1]. Then 
for all g E C\j we have 

Pci{9)= PcA9) = \\9\\op- (3.35) 

This demonstrates that £^ is symmetric, since we have 

Pci,{9) = p{9) = \\9\\op (3.36) 

for all g G £^, and thus of course (3.36) holds for self-adjoint elements as well, which 
implies satisfaction of Hulanicki's lemma for C\,. We can also conclude, using distri- 
butions closely related to pseudoinverses in a technique which will be described later 
in this section, that C\, is inverse closed in B{LF'{W^)) for the above weights. It is 
apparently an open question whether or not C^ is symmetric (see Section 8.2). 

We restate a lemma from [23]: 

Lemma 3.3.4 (Grochenig, Leinert, Pytlik) For any unbounded admissible weight 
function v there exists a sequence of admissible weights Vn, nG N, with the follow- 
ing properties: 

(a) Vn+l <Vn<v, \fn 

(b) 3 c„ > such that v < CnVn 



-^W^ni 
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(c) lim^^oo'^n = 1 uniformly on compact sets ofW^^ and 

(d) 3 I3n such that Vn{z) = e~"'v{z) for \\z\\ > /9„. 

3.3.3 Cl is inverse-closed in B{L'^{R'^)) 

Lemma 3.3.5 Suppose v{x) > C(l + \x\Y with 6 > 0. Under the hypotheses of 
Lemma 3.3.4 the following identities hold \/ f = f* E Cl : 



(a) lim||/|Uj^ = ||/|Ux, and (3.37) 

(b) pc.if)=Pc<f) = \\f\\op (3.38) 
Proof, (a) Let e > 0. Then if / = /* = / + /le G Cl, we have 



ess sup / \f{x,y)\vn{x-y)dx. 



Part (d) of Lemma 3.3.4 states Vn{z) = e "-v^z) for ||2;|| > /9„, so 3 no = no(e) G N 
such that 



ess sup / \f{x,y)\vno{x-y)dx<e """WfWci <e. 

yGRd. J 

xeR'^:\\x-y\\>f3„o 

Invoking this and monotonicity Vn+i < Vn < v, we get for all n > Uq, 



ess sup / \f{x,y)\vn{x—y)dx<e. 

By property (c) of the f„, we have that f „ -^ 1 uniformly if ||2;|| < /3„q, so for 
n > ui = ni(e), 

ess sup / \f{x,y)\vn{x-y)dx < (1 + e) ess sup / \f{x,y)\. 

xeK'';||a;-?/||</3„Q 

These last two inequalities yield, for n > max{no,ni}, 

||/|Ui^<e + (l + 6)||/|Ui. (3.39) 

We take the limit of (3.39) to get lim ||/|Uj < ||/|Ui. Since ||/|Ui =\\f\\ci +|/i| and 
WfWci = WfWci + lAt|, we obtain lim ||/||£i < ||/||£i. It is clear that since Vn, w„ > 1, 
we have lim ||/||£i > \\f\\c^, so we conclude lim ||/||£i = ||/||£i. 

n—>oo " n— >oo " 

(b) Since v and f„ are equivalent and p{g) < \\g\\ for any Banach space norm, we have 

PciUf = Pciif) = PciSf) < Wf'Wci^ VA;,n,GN. 
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This and part (a) yield 

PciUf< lim||/'^|Ui^ = ||/|U. VA:,GN. 

Taking k^^ roots gives, by continuity of x'^, 

PciU)< lim||/1iy.' =PcAf)- (3.40) 

n— >oo -"n 

Now v{x) > C(l + |a;|)'^ = ts{x) and 6 > ^ C^ C C^. By Lemma 3.3.1, we have the 
reverse inequahty of (3.40), namely Pc^{f) < Pciif)- Thus Pc^{f) = Pciif)- Barnes' 
Lemma 3.3.3 implies pc^{f) = \\f\\op- We conclude pciif) = WfWop- □ 

This is an important lemma, since it affords us the following theorem. 

Theorem 3.3.6 Assume thatv is an admissible weight which satisfies the weak growth 
condition 

v{x)>C{l + \\x\\Y for some 5 >0. (3.41) 

Then V / G Cl, o'£i(/) = o'(/), thus Cl is a symmetric Banach algebra and inverse 
closed m B{L'^(R'^))^. 

Proof. By (3.38), Cl satisfies the hypotheses of Hulanicki's Lemma 3.3.2, and we 
obtain crciif) = ^(Z)) ^ f = f* ^ '^l- Thus is Cl a symmetric Banach algebra. Also, 
if / = /* G £^ is invertible on L^(]R'^), then / is invertible as an element of Cl- To 
apply this conclusion to elements of Cl which are not self-adjoint (i.e., to conclude 
that Cl is inverse closed in B{L^{R'^))), let / G i:^ be invertible on B{L^{R'^)) but 
not necessarily self-adjoint. Since / is invertible on B{L'^{M.'^)), we know /*/, //* are 
invertible on i3(L^(]R'^)) as well. Indeed, if / is invertible 3 g 3 fg = gf = e = 
S{x — y). Now 5{x — y) is self-adjoint, and by (2.25) we have 

g*r = if 9)* = S{x -y) = {gfr = fg* 
so /* has inverse g*. Then by (2.23) 

ifrwg) = u\rf)9) 

= f{S{x-y)g) 
= f9 = S{x-y); 

the inverse of /*/ is computed similarly. We conclude that f*{f f*)~^ and {f*f)^^f* 
are, respectively, right- and left-inverses of / in Cl- Thus even if / is not self-adjoint, 
if / is invertible on B{L'^(M.'^)), then / is invertible on Cl, and thus Cl is inverse closed 
in i3(L^(]R'^)). We conclude that if / G Cl, then the right- and left-inverses of / share 
the same decay properties as /. □ 
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3.3.4 Cy is inverse-closed in B{L'^{W^)), for s > d 

We next consider Eu^s and C^. The main result: 

Theorem 3.3.7 Assume that u is an admissible weight, s > 0, and v = UTg. Then 

PBUf) = \\f\\op y f = reBu,s. (3.42) 

Thus Bu,s is a symmetric Banach algebra, and inverse closed in i3(L^(]R^)). 

Proof. Writing v = {uts'jts., we can assume without loss of generality that u satisfies 

condition (3.41), i.e., £i C -C^ /2- Let f = f + fie E Bu^s C £„. Then (3.28) yields 

ii/'ik„<2K|iriUiiiriu„, vn>i. 

which then implies 

WflKs = ll/%:.+^"IUl+j^l'" 

< 2«:||7^|Ux||7^|U„+2«:||7^||^,+2«:|/ip- 

< 2K{\\nci^ + \f,n{\\nci + \\nc. + \f^n 

= 2K\\nci\\nB.^^ (3.43) 

The spectral radius formula is PA{g) = li^i US'"!!//"- So we have 

n—>oo 

PB.M) = i^jni'C 



= \IpciU')\IpbuM')^ 

which implies 

pBu,s{f) <PClU')- 

Since £^ C jCIj^, we can apply (3.38) of Lemma 3.3.5 and get 

PBUf) < PClif) = WfWop. 

It is clear that 

pBuAf) >pciU) 
since Bu,s C Cl^. Thus we have shown 

PBuM) = Pciif)- 

We conclude that o'b^^(/) = cr{f) and Bu,s is symmetric and inverse closed in 
i3(L2(M^)). "'' ' □ 
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Finally ii s > d, the algebras Bu,s and C^ coincide, and we have the following: 

Corollary 3.3.8 Suppose thatu is an admissible weight, s > d, and v = uTs. Suppose 
also that f is invertible on L^(]R^) and satisfies the ojf- diagonal decay condition 

ess sup\f{x,y)\v{x — y) < oo. (3.44) 

Then g = f^^ satisfies the same condition, i.e., 

ess sup\g{x,y)\v{x—y) < oo. (3.45) 



Proof. Suppose first that /* = /. Since s > d, Bu^s = '^v and Theorem 3.3.7 says 
(7Bu,s{f) = ^{f)- Thus since ^ o'(/), we know that ^ <^t3u,sif)- ^ow since / is 
invertible in £„, the inverse of / must share the defining property of C^, namely 
(3.44). To extend the result to all f E Cy, we apply the same argument as in the 
proof of Theorem 3.3.6. □ 

We also note that if / G £^ or / G £^, then in the case that / is a positive operator, 
all powers /", a G M are in Cl or C^, respectively, by the Riesz functional calculus 
(see e.g., [47]). 



Template of an example 

As an illustration of how a specific example of an operator in Cl or C^ might be 
discovered, we investigate the convolution operator generated from f{y) = S{y) + 
24 e~''^^ Note that we have / = ! + /. Let W = W(M.'^) denote the Wiener algebra, 
for which g eW ^ g,g e L^{W'-), and fix h e W. Then 



^iii = / (i+/Mr 



/, 



hiy)e-''''--ydy 



e'^^-^rfcu (3.46) 



%) / il + fiu)r'e'''''"<''-y^du;dy (3.47) 



Hy)J'-' I ~ J (x - y) dy (3.48) 

h{y)J^-' I 1 - -S^ ] {x - y) dy (3.49) 



M.)^-'(i-^) 



[x - y) dy (3.50) 
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h{y) 



5{x -y) -J^ 



-1 



(x-y) 



dy (3.51) 



the non-unit part of the kernel is f~^{x,y) = —T 



-1 



^ — I \x — y\. Since for all 



o; G M , we have 



< 



1 



< e~ 



we conclude that T 



-1 



i^'^'^^l 



2 - e'^^' + 1 

shares the same decay as 23 e" 



(3.52) 



T-H2ie 



Of course e '"^ is a member of neither C\ nor £^, but finding an explicit example in 
C\ or £„ has thus far proven elusive. 
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Chapter 4 

Pseudodifferential Operators and 
Linear Time- Varying Systems 



The results in Chapter 3 should find application in the area of wireless communica- 
tions channels. Let s{t) be a wireless signal transmitted through a channel H, and 
let s{t — t) denote its r-delayed version. The received signal is 

r(t) = (Hs)(t) + e(t) = / h{t, T)s{t - T)dT + e(t), (4.1) 

where /i(t, r) is the impulse response of the time-varying channel H [3] and e(t) is 
zero-mean stationary white noise (4.1 will be explained fully in Chapter 5). A typical 
channel may have a line-of sight component, and cause many refiections from objects 
in the vicinity of the transmitter and/or receiver (as well as Doppler effects from 
the relative motion between the transmitters and receivers). Thus H will yield an h 
which consists of the sum of 5{t — r) and some other function of t and r possessing 
off-diagonal decay, quite possibly granting h membership in C\j or C^. The wireless 
channel applications of the work in Chapter 3 will be the subject of future investiga- 
tion. 

There is an interesting connection between pseudodifferential operators and linear 
time- varying systems, of which a wireless channel is one example (see [57]), that will 
be instructive to describe as a point of departure for our excursion into signal design 
for wireless communications. 



4.1 Pseudodifferential operators 

We will assume that all functions described in this section are in 5(]R'^) so that the 
integrals described herein behave; i.e., we interpret all expressions in the distributional 
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sense. In the study of partial differential equations one considers equations of the form 

Af{x) = J2 <yc.{x)D-f{x) = g{x) (4.2) 

|a|<Af 

where A^ is the order of A, a is a multi-index, and {(J a] is the set of non-constant 
coefficients. We have 

{Dy){uj) = {2mujrf{uj) (4.3) 

(see section 2.5). Using (2.13) and (4.3), one can write (4.2) as 






|/(cu)e'"'"-"rfcu (4.4) 

"■\a\<N 

a{x,u)f{u)e^^'''-'' du (4.5) 

where a is the symbol of A. One can admit more general symbols, yielding the defi- 
nition of a pseudodifferential operator. 

Definition 4.1.1 Let o" be a measurable function or a tempered distribution on M?'^. 
Then define 

KJXx):= [ a(a;,cu)/(cu)e2— rfcu. (4.6) 

The operator K^r is called the pseudodifferential operator with symbol cr. 

With suitable restrictions on a, K'^ becomes a multiplication operator, or a Fourier 
multiplier, or a differential operator with constant coefficients. Additionally, we can 
write K'^ in several ways. First, use (2.13) in (4.6) to yield 

KJXx) = f ( [ a{x,uj)e-'-'^y-^>-duj)f{y)dy (4.7) 

K{x,y)f{y)dy (4.8) 

where we deonte by k^r the kernel of the integral operator Krj. Note then that the re- 
sults of Chapter 3 apply to pseudodifferential operators. Denote by La the coordinate 
transformation LaF{x,y) = F{x,y — x) and let JF2 be the partial Fourier transform 
in the second variable. Then k^r can be written as 

k{x, y) = J^2cr{x, y - x) = LaJ^2<y{x, y). (4.9) 

We can also interpret (4.7) as the convolution of f{y) with J^2^a{x,y) from 

K,f\x) = [ .F2- V(x, X - y)f\y) dy. (4.10) 
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4.2 Linear time-varying systems 

A linear system transforms an input signal / into an output signal Af for which 
^(ci/i + C2/2) = CiA/i + C2A/2 for all inputs fj and constants Cj. For a time-invariant 
system, we have A{Txf) = T^Af V a; G M"'. Any linear operator satisfying this com- 
mutation relation can be written as a convolution operator 

Af = f*h where h = A6. (4.11) 

The frequency response of the system A is 

A(e2™-*)(a;) = / /i(t)e=^™-(^-*) rft = /i(cj)e2™-^ (4.12) 



and we conclude that the complex exponentials e'^'^^'^'^ are the eigenf unctions of A 
with eigenvalues h{uj). 

A linear system A may change with time, in which case A is called a time-varying 
system. In this case the convolver h depends on time, and we have 

Af{x)= [ h{x,x-y)f{y)dy. (4.13) 

By (4.10) we can interpret the linear time- varying system A as a pseudodifferential 
operator K^r with the symbol a = J-'2h. 

In general one has K~^ 7^ -^i/o-, but for many pseudodifferential operators, in partic- 
ular underspread operators, where the the channel does not spread the signal "too 
much" in the time-frequency plane, (see Chapter 5), one has K~^ ^ -ft'i/o- (see [32]), 
and this is the case for a typical wireless communications channel. Underspread 
operators have been studied, but these have almost always been Hilbert-Schmidt 
operators, which are not invertible. This is too restrictive for the case of wireless 
communications, since it excludes the case where the channel operator is just a small 
perturbation of the identity. 

A wireless communications channel is a linear time-varying system, and can thus 
be interpreted as a pseudodifferential operator. Denote by Hg. the channel whose 
characteristics are described by the symbol a. Strohmer proved in [57] that if o" is a 
wireless channel operator symbol, then a G M^'^, the space of symbols a for which 

snpV^aiX, n)m{n) dn < 00 (4.14) 

X 

where X, fi G M^*^, m is admissible in the sense of Section 3.1, and V denotes the 
2(i-dimensional STFT. 
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In the case that the channel is underspread and the pulses used for encoding the data 
have good time-frequency localization or TFL (explained more rigorously in Chapter 
5; for now we will simply say that a pulse has good TFL if it possesses good decay 
in both the time domain and frequency domain), we have H~^ ^ Hi/a in the weak 
sense, meaning that we can trade invertibility for approximate invertibility and sig- 
nificantly less effort. 

In the sequel we shift the focus from the operator to the functions which make up the 
components of our signal. To ensure that the spreading induced by a given channel 
operator is minimized, we consider, and make some demands of, the TFL of the pulses 
on which the data is encoded and sent through the channel. 

In this regard, we have developed several improvements to a communications scheme 
called orthogonal frequency-division multiplexing, or OFDM. This system not only 
shows great promise as a candidate for future wireless communications systems, but 
possesses a beautiful mathematical foundation as well. The remainder of this thesis 
will be concerned with our improvements to the algorithm currently used in practical 
applications of OFDM. 
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Chapter 5 



OFDM 



5.1 Introduction 

OFDM is an efficient technology for wireless data transmission [61], and is currently 
used in the European digital audio broadcasting standard [10], in digital terrestrial 
TV broadcasting, broadband indoor wireless systems [38], and an increasing number 
of wireless internet service providers (WISPs) in the U.S. [49]. The advantages of 
OFDM include robustness against multipath fading and high spectral efficiency. 

Transmission over wireless channels can be subject to multipath propagation (i.e., 
that the received signal is the sum of perhaps a line-of-sight component and several 
reflected delayed components), leading to time dispersion which causes intersymbol 
interference (ISI), and the Doppler effect (frequency compression or rarefaction due 
to the constancy of the speed of light and the relative motion between the transmitter 
and receiver) which causes dispersion in the frequency domain leading to interchannel 
interference (ICI)[41]. The distortion resulting from time-frequency dispersive chan- 
nels depends crucially on the TFL of the pulse. For instance OFDM systems based on 
rectangular-type pulses with guard-time interval or cyclic prefix can prevent ISI, but 
can not effectively combat ICI, since they decay too slowly in the frequency domain 
(see Figure 5.2). 

The design of time-frequency well-localized pulse-shaping OFDM filters is therefore 
an active area of research [36, 58, 24, 32, 8]. Some of the proposed pulse shape design 
methods involve the solution of highly non-linear optimization problems [58, 24], 
other approaches lead to pulse shapes that can be easily computed, but little or no 
theoretical justification of the optimality of the so-obtained pulse shapes has been 
given [36, 8, 32]. 
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5.2 OFDM basics 

The principal idea of OFDM is to divide the available spectrum into several subchan- 
nels (subcarriers) so that over each subcarrier the data rate is slower, the symbol 
period is longer and thus the system is more tolerant of time dispersion. 

The transmission functions of an OFDM system consist of translations and modula- 
tions of a single pulse shape ip, i.e., 



iJkiit) = ij{t - kT)e 



2mlFt 



(5.1) 



for some T,F > 0,k,& Z and i = 0, . . . , N — 1 or i ^ Z. We denote such a function 
system (whether orthogonal or not) generated via (5.1) by the triple (■?/', T,F). We 
call T the symbol period and F the carrier separation. Let W be the given bandwidth, 
then Nc = W/F is the number of subcarriers. Throughout our analysis it is convenient 
to assume A^ = oo, so the frequency index i in (5.1) runs through the integers instead 
of from to A^ — 1. But the results presented in this paper also hold for finitely many 
subcarriers, either verbatim or with a negligible approximation error. A baseband 
OFDM system is schematically represented in Figure 5.1. 



Transmitter 



Channel 



Receiver 



(k+I)T 




(Hs)(t) 




Figure 5.1: A general baseband OFDM system configuration 

A necessary condition for perfect reconstruction of the data is that the functions i/jk/ 
are linearly independent. This implies that TF > 1 which follows easily from the 
connection between OFDM systems and Gabor frames, e.g., see Section 6.5.1 of this 
paper or [32]. We note that orthogonality of the functions ipk/ is not a requirement 
for perfect reconstruction, but minimizes the error caused by AWGN. 



Let Ck/, /c, £ G Z, be the data symbols to be transmitted. The OFDM signal is then 
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given by 

s{t) = y^ Ck,eiJk,e- (5.2) 

k,e 

The received signal is 

r(t) = (Hs)(t) + e(t) = / h{t,T)s{t-T)dT + e{t), (5.3) 

where h{t, r) is the impulse response of the time-varying channel H [3] and e(t) is 
zero-mean stationary^ white noise. The data symbol Ck/ is transmitted at the {k, £)*'' 
point of the rectangular time-frequency lattice {kT,iF),k,i G Z. Assuming that 
the ipk/ are orthonormal the transmitted data are recovered by computing the inner 
product of the received signal r with the functions ipk/, i.e., Ck/ = {r,ipk,e)- By (5.2) 
and (5.3) we have 

Ck,i = {Us, i)k/) + (e, i)k/) = ^ Ck'/' {Uipk'/', ^k,e) + (e, ^k/), (5.4) 

k',i' 

where the interchange of summation and integration is justified under mild smooth- 
ness and decay conditions on ip and H. In the ideal (but unrealistic) case the wireless 
channel does not introduce any distortion. Thus the orthogonality between the func- 
tions ipk/ is preserved, and if e = the data Ck/ can be recovered exactly. However in 
practice wireless channels are time-dispersive and frequency- dispersive. Hence when 
ipk/ passes through the time-varying channel H it is spread out in time and fre- 
quency, causing ISI and/or ICI in the transmitted data. It is customary (cf. e.g. 
[36, 41, 24, 32]) to describe the channel by the maximal excess delay Tq and the m,ax- 
imal Doppler spread z/q of the channel scattering function. It is clear from (5.4) that 
the interference of the data symbol Ck/ with neighbor symbols will be smaller with 
better TFL of ipk/ and with larger distance between Ck/ and adjacent data symbols 
in the time-frequency domain. 

As shown in [24] for a large number of subcarriers the spectral efficiency p of the 
OFDM system in terms of data symbols per second per Hertz^ is approximately 
given by p = {TF)~^. Since TF > 1 we have of course that the maximal spectral 
efficiency is Pmax = 1- Thus any increased distance between adjacent symbols in the 
time-frequency domain costs the designer in terms of spectral efficiency. 



^A stationary stochastic process is defined by the statistical invariancc of its properties under a shift of the tiine- 
Lgin. 

^In the remainder of this paper we will use p to denote specti 
will be no confusion with the spectral radius p(f) of Chapter 3. 



origin. 

^In the remainder of this paper we will use p to denote spectral efficiency, since it is standard notation, and there 
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5.3 Dispersion in the wireless channel 

5.3.1 The Doppler effect and ICI 

The formula for the Doppler shift is given by 



fd = -cos(6') 



(5.5) 



where fd is the change in frequency of the signal due to v, the relative velocity of 
the transmitter and receiver and A, the wavelength of the signal (for simplicity, we 
assume a pure frequency), and 6 is the instantaneous angle of approach. For a speed 
of 1100 mph (^492 m/s, the approximate relative velocity of two approaching airlin- 
ers), on direct approach {6 = 0), for a 5 GHz signal (A = 6 cm) the Doppler shift 
is then 8.196 kHz. But for IEEE 802.11a, the subcarrier spacing is 312.5 kHz. So 
how can such a small relative shift in frequency lead to significant ICI at the receiver? 




Figure 5.2: Frequency domain representation of five standard (square-wave-carried) OFDM tones. 
Note that the peak of each tone is at a zero of every other tone. 



Figure 5.2 illustrates the frequency domain representation of the modulated square- 
wave pulses used in standard OFDM. Note that the peak amplitude of each frequency 
is located at the zeros of all other frequencies. This is what distinguishes OFDM 
from FDM, a system wherein the designer precludes interference only by maintaining 
a large guard interval between frequencies, which can reduce spectral efficiency dra- 
matically. OFDM requires no frequency guard interval, since ideally the subcarriers 
do not interfere with each other. But as one might suspect, this system is quite sen- 
sitive to frequency perturbations. When the receiver samples at the center frequency 
of each subcarrier, it gathers only energy from the desired signal, provided that there 
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is no frequency dispersion. This is referred to as zero-delay orthogonality. It is clear 
from Figure 5.2 that it is not necessary for the Doppler shift to be anywhere near half 
of the subcarrier separation to cause significant interference. One cannot speculate 
as to the relative velocities which communications systems might encounter in the 
future. 

5.3.2 Multipath delay and ISI 

A typical received signal consists of (perhaps) a direct line-of-sight signal from trans- 
mitter to receiver, plus reflected signals arriving at later times due to the necessarily 
greater path lengths travelled, spreading the signal energy in time. Figure 5.3 shows 
this effect. 

Transmitted Signal 




Reflected Signal 







Received Signal 




Figure 5.3: Multipath delay spread affects a signal in a reflective environment. 



5.4 Time-frequency localization 



Ideally we would like to construct an OFDM system (■?/', T,F) that satisfies the fol- 
lowing three requirements simultaneously: 

(i) the functions ipk/ are mutually orthogonal. 
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(ii) ip is well-localized in time and frequency domain, 

(iii) TF = 1, i.e., maximal spectral efficiency. 

As is well-known from Gabor theory these three conditions cannot be satisfied si- 
multaneously due to the Balian-Low Theorem [14] the typical statement of which 
precludes, in essence, simultaneous satisfaction of (ii) and (iii). It should be noted, 
however, that this holds true, even if we relax condition (i) and consider biorthogo- 
nality instead (which however would increase the sensitivity of the system to AWGN). 

Hence we relax condition (iii) and allow TF > 1. This loss in spectral efficiency 
is usually an acceptable price to pay to mitigate ISI/ICI (provided TF is not too 
large). As pointed out above in a dispersive channel the transmitted signal is spread 
in time and in frequency, causing ISI and ICI. This interference can be reduced if 
the pulse shape possesses good TFL. The classical way to measure the TFL involves 
the Heisenberg Uncertainty Principle. More precisely, if / G 2y^(]R) and r, cu G M are 
arbitrary, 

(t - rr\m\'d?j ' ^(z. - c.)V>)Prf^) ' > ^. (5.6) 

Equality holds if and only if / is of the form f{t) = c^'^^^'-'it-'^) e-^^^i^-^) for t,uj E 
M, c 7^ 0, and a > 0, i.e., if / is a (possibly translated and/or modulated and/or 
dilated) Gaussian. 

In (5.6) time and frequency are treated separately. It seems more appropriate in our 
context, in particular with respect to time- frequency dispersive channels, to measure 
the TFL by considering time and frequency simultaneously. Moreover in light of 
(5.4) we are naturally interested in measuring the cross-interference between differ- 
ent transmitter functions ipk/ and ipk',i'- Thus we replace the pair (/, /) by a joint 
time-frequency representation of two functions / and h, such as the cross-ambiguity 
function A{f, h). 

In light of the Heisenberg Uncertainty Principle it is not surprising that A{f, h) cannot 
be arbitrarily well concentrated. This can be expressed for instance by the fact that 
A{f,h) cannot be compactly supported, see [28]. An alternative formulation is the 
following Heisenb erg-type inequality [11, 9]. 

Theorem 5.4.1 For every f,hE X^(]R) there holds 

it-Tf\A{f,h)it,u)\^dtdu- f f{u-ujf\A{f,h){t,u)\^dtdu> 




47r2 ■ 
(5.7) 



CHAPTER 5. OFDM 43 

Equality holds if and only if 

f{t) = cie2"*'^*e-™(*-^)', 
h{t) =c2e'"*("+'^)*e-™(*-"-'*)', (5.8) 

with a > 0, Ci, C2 7^ 0, and P,s,r E M. 

Thus again, translated, modulated, and dilated versions of the Gaussian lead to op- 
timal TFL. But Gaussians cannot be orthogonal under translation, so we appeal to a 
technique first discovered in the finite-dimensional discrete case by the chemist Lowdin 
and here termed Lowdin orthogonalization. By applying an infinite-dimensional con- 
tinuous version of Lowdin's technique to form orthogonal pulses from the Gaussian, 
we retain as much of the Gaussian's optimal TFL as is theoretically possible. 
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Chapter 6 

Frames, Gabor frames, and Lowdin 
orthogonalization 



To construct an OFDM communications system which is robust against both ISI 
and ICI, one must begin with a function h possessing good decay in both time and 
frequency, and generate a set {h^/} of functions by integral time- and frequency- shifts 
of g. This setup is almost identical to a general Gabor system (defined in Section 
6.3.2), with the exception of the phase factor present in a Gabor system due to the 
application of modulation first, and then translation to the Gabor system generator, 
but with the additional constraint that every function in the sequence {hk/} be 
orthogonal to all other members of {hk/}- Since any system {gu/} generated from 
the Gaussian g obviously does not yield an orthogonal set of transmission functions 
defined in this manner, and nor do any dilations of the Gaussian^, we use a standard 
method to "orthogonalize" the set {gu/}- In connection with pulse shape design it 
has been used implicitly in [36] (see also [6]), and explicitly in [7, 8, 54]. We describe 
this algorithm first in finite dimensions as an illustration of its elegance. 



6.1 Discrete Lowdin orthogonalization in finite dimensions 

Suppose $ G C™"^", with m > n and rank($) = n. Consider a nonsingular linear 
transformation A G C"^" which, when applied to $, yields an orthonormal basis 
^ = •I'A of C". We have 

/„^„ = ^*^ = (<|.A)*(<I>A) = A*$*$A, (6.1) 

and we define A G C"^" as the Gram matrix $*$. By the properties of (■,■), A 
is Hermitian, and by Proposition 5.1.5 in [22], A is a positive invertible operator. 



-"^ Since for all cr > , g" {t) = (2o-)^'''e '^'^' > for all t, gcr{t) can never by orthogonal to any of its shifted copies, 
a similar argument holds of course for modulations of g„- . 
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Consequently, 3 U unitary which diagonahzes A := U*AU, and the elements of A are 
positive. We can define an arbitrary function / on A via 

/(A) = U*f{A)U, (6.2) 

provided that the Taylor series for /(A) converges, since (6.2) is always true for poly- 
nomial /. In particular, for f{x) = a;"^/^, we note that A is nonsingular, so the 
Taylor series for A^^/^ converges just as the series for f(x) = x'^^"^ converges for 
X > 0. I.e., if a; G (0, oo], then since the mapping f : x ^ f{x) is continuous and 
infinitely differentiable on (0, oo], the Taylor series for / converges at x. Thus we have 
A-i/2 ^ f/*A-i/2f/^ and A'^/^ jg Hermitian as well. 

Define B = A^/^A, so that A = A'^/^B and A* = B*{A-^/^)* = B*A-'^/^. Thus 

A*AA = B*A-^/^AA-^/'^B = B*B (6.3) 

and we see that the equation A* A A = I has a solution if B is an arbitrary unitary 
matrix. The choice oi B = I leads to Lowdin or symmetric orthogonalization (or 
orthonormalization): ^ = $A^^/^. 

An important property of Lowdin orthogonalization is that of all possible orthogo- 
nalizations it disturbs the original basis the least under the Frobenius norm, defined 

by 






|2 



Theorem 6.1.1 Let ^ be a linearly independent set of n vectors in C™". Considering 

$ as an mx n matrix of complex numbers, we have the singular value decomposition 

$ = UDV*, where U G C™^", V G C"^", U and V are unitary, and D G C"^" 

contains the singular values of ^. Then if Q is unitary, ||$ — Q\\f is minimized by 
Q = P= $($*$)-i/2 = f/y* 

Proof. The strategy is to show first that P = UV* and thus that P is unitary, and 
then to show that P minimizes the above expression. 

$ = UDV* =^ $($*$)-i/2 = UDV*{VDU*UDV*)-^/^ (6.4) 

= UDV*{VD^V*y^/^ 

= UDV*{VDV*VDV*y^/^ 

= UDV*{VDV*)-^ 

= UDV*VD-^V* 

= UV* 
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Thus P = UV* and P is unitary. To show that Q = P minimizes ||$ — QHf among 
unitary Q, note that 

11$ - QWf = \\UDV* - QWf = \\D - U*QV\\f, 
which is minimized if U*QV = I, or Q = UV* = P. □ 

Having treated the discrete finite case, we turn to the continuous infinite case, begin- 
ning with an introuction to frames. 

6.2 Frames 

There are many ways to reconstruct L^(]R'^)-functions using sets of functions in 
L^(]R^). We can use wavelets, or even more exotic sets of functions. Orthonormal 
bases are usually ideal in some sense, but often enough we have no stake in unique- 
ness of coefficients. This particular indifference leads to the use of what are called 
frames. A sequence {ejjjgj, over some index set J, in a Hilbert space 7i, is a frame 
for 7i if and only if 

3A,BeR+3yfen AWff < J2\{f,e,)\' < B\\ff. (6.5) 

For any subset {ej}j(zj C 7i, we have the analysis and synthesis operators 

Cf = {{f,e,)\3eJ} 



DC = Y. 



CjCj. 



It can be easily shown that D = C* . The action of the frame operator S = C*C = 
DD* on f eHis 

Sf = Y.{f,e,)e,. (6.6) 

The condition A > ensures that S is injective; the condition B < oo ensures that 
S is bounded and thus continuous (see any elementary text discussing operators on 
Banach spaces). Since {Sf, f) = {C*Cf, f) = {Cf, Cf) > 0, S* is a positive operator. 
li A = B then {cj} is a tight frame. Tight frames are important to us because they 
can lead to orthonormal systems, which we will see later. 

Lemma 6.2.1 If {cj \j E J} generates a frame, then the set {S^^Cj \j E J} gener- 
ates a tight frame. 
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Proof. (From [22]) 

/ = S-^^SS-'^f = 5^(/, S-h,)S-h, (6.7) 

^ {fJ) = J2\{f,S-h,)\'. (6.8) 



We will refer to frames {dj | j G J} and {cj | j G J} for which 

/ = $^(/,e,)rf, (6.9) 

iGJ 

iG"' 
as dual frames with respect to the identity, or simply dual frames. If 

S'-f = 5^(/,e,)rf, (6.11) 

iGJ 

= E(/'^^>^- (6.12) 

iGJ 

then we refer to {dj} and {e^} as dual frames with respect to 5°. 

6.2.1 Powers of the frame operator 

Assume that {cjjjgj generates a frame. Consider the nonorthogonal expansions 

/ = J2^f^S-h,)S~'^ej (6.13) 

iGJ 

= 5^(/,e,)5-S. (6.14) 

JGJ 

= E(/,5-S)e,. (6.15) 

JGJ 

We can use these to derive several results about the frame operator S . 

• For which p does S^Cj generate a frame? It follows from Appendix A. 7 in [22] 
that VpGM, S Cj generates a frame. 
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What are the frame elements which generate the frame operator 3 7 If we 



operate on / in (6.13), we have 



q—r 



57 = s {s-" /) 

= J2(s''f,s-h,)s'-\s'hj) 

= Y.^f , S-K,)S'-'-K, (6.14) 

where we used the fact that powers of S are self- adjoint, and the expression 
(6.13) (as well as (6.14) and (6.15)) converge unconditionally. We conclude that 

the S^~ej generate the frame for S {r = q/2), and that the S "^ Cj and S Cj 
generate frames which are dual with respect to S . Equivalently, if the frame 
elements are 5**6^, then the frame operator is 5'*^^*"''^^/^. 

• What is the dual frame for {S Cj^j^jl We have 

/ = S-\s'f) (6.15) 

jeJ 
= Y.{f,S''e,)S-'^\, (6.15) 

which shows that the dual of S* is S* 

6.2.2 Riesz bases 

Let J be some index set. 

Definition 6.2.2 A set {ej}j^j is a Riesz basis for 7i -^ 

3A',B' e M+ such that A'||c||^ < ||5^9ej||^ < 5'||c||^ (6.15) 

for all finite sequences c. 

What is the essential difference between a Riesz basis and a frame? A Riesz basis 
yields unique coefficients while a frame does not; alternatively a Riesz basis is a 
spanning set of minimal cardinality while a frame is redundant. It is instructive 
to present an example distinguishing frames and Riesz bases. Consider the finite 
Hilbert space M^ and the frame (but not Riesz basis) Gi = {(1,0), (0, 1), (1, 1)}. Gi 
is not a Riesz basis since we cannot find an A' > for which (6.2.2) holds (example: 
Ci = C2 = 1, C3 = — 1 yields an inner sum of 0). But Gi is a frame since the inner sum 
in the inequality (6.5) is positive V/ G M^, and obviously finite. Clearly any Riesz 
basis is a frame since we have uniqueness of coefficients and can take Cj = (/, Cj). 
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6.3 Lattices and Gabor frames 

The main distinction between a frame and a Gabor frame is that a Gabor frame is 
necessarily defined on a symplectic lattice in M^'^. 

6.3.1 Lattices, symplectic lattices, and adjoint lattices 

We briefly review some basic facts about symplectic lattices and lattices in general, 
see [12] for details. The generator matrix L of a lattice A contains the coordinates of 
the column vectors that describe the fundamental domain of a lattice; i.e., A = L{7/). 
The density ^(A) of a lattice A is given by ,^(A) = l/det(L). Write A = (x, a;) for 
x, tu G W^. The symplectic form is [Ai,A2] = X2 ■ oJi — xi ■ uj^- A lattice A with 
generator matrix L which preserves the symplectic form, i.e., [LAi,LA2] = [Ai, A2], is 
called symplectic, and is an element of the symplectic group Sp{d). If A G Sp{d), then 
det(L) = ±1, so Sp{l) =SL{2,R). 

Definition 6.3.1 Any lattice A whose generator matrix L = cA , A E Sp{d) , c eM., 
is called a symplectic lattice. 

Definition 6.3.2 The adjoint lattice of a symplectic lattice is A° = ,^(A)~dA, where 
L is the generator matrix for the original lattice A. 



A separable or rectangular lattice has generator of the form L 



a,f3 e Z'^. 



'a 0" 

p_ 

In the sequel, all lattices under consideration will be symplectic, unless otherwise 
indicated. 



6.3.2 Gabor frames 

Let Tx denote the translation operator whose action on / is given by Txf{t) = f{t—x), 
and denote by M^ the modulation operator, the action of which on / is given by 
Muifit) = /(t)e^'^*'^'*. We then have the commutation relations 

T,M^ = e-^^'^-'M^T,, (6.15) 

and it is quite easy to demonstrate that 

{TJ) = M_J (6.15) 

and 

{MJ)=tJ (6.15) 

Definition 6.3.3 Fix both nonzero g, G LF'{R.'^) and a lattice A, with generator ma- 
trix L and det(L) > 0. Write z{\) = T^M^^ , A = {x,u) G A). The set of time- 
frequency shifts 

{g,A) = {z{X)g\XeA} (6.15) 
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is called a Gabor system, or Gabor frame if a frame. 

a 



p 



a,p e 



If 7 = (yf we say that {g, A) is a Gabor frame. If A is separable with L 

7j'^, then we write 

{g, a, /3) = {gke := T.^M^.^? | A;, £ G Z'^}. (6.15) 

Note the similarity of an OFDM system to a Gabor frame; the structure is very 
similar. A Gabor system includes a phase factor which occurs since modulation is 
performed first, then translation, while the order of operations is reversed for OFDM, 
so for OFDM there is no phase factor. If the product x ■ cj G Z, the phase factor 
disappears. Aside from the phase factor mentioned at the beginning of this Chapter, 
there is additional difference between Gabor frames and OFDM systems, namely that 
they are dual in a certain sense. Analysis and synthesis by a Gabor frame is 

/ G L'{R) ^ {dki} G f{Z') ^ / G L'm (6.15) 

where the frame must be complete (i.e., the closed linear span of the frame is the 
entire space) and may be linearly independent for perfect reconstruction of /. OFDM 
eoncoding, transmission, and data recapture is 

{Cki} G f{Z') ^ s G L\R) ^ {dki} e f{Z') (6.15) 

where the ipki, of which s is a linear combination, must be linearly independent and 
may be complete for perfect reconstruction of {cm}- 

Proposition 9.4.4 in [22] states that a Gabor system on a symplectic lattice is the 
image of a Gabor system on a separable lattice under a unitary operator. 



Definition 6.3.4 The Gram, matrix Rg,A associated with {g, A) is defined via 

{Rg,A}x,\' ■■= {9x',9x)- (6.15) 

If A is rectangular, we write Rg,T,F- Rg,k may be invertible, as is the case when {gx} 
is a Riesz basis, or singular, as when {g\\ is a redundant (possibly tight) frame. 

The action of the Gabor frame operator S is: 

Sf = Y,U\^W9)^{^)9 (6.15) 

AeA 
or 

Sf = ^ {f,9ke)9ke (6.15) 

if A is separable. If we wish to emphasize the dependence of S on the function g and 
the lattice A, we write Sg^g^A- 
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Definition 6.3.5 Let (7,7 G L^(]R'^), and suppose that every / G L^(]R'^) can be 
expanded as 

/ = J2(f,z{X)g)z{Xh (6.16) 

AeA 

= 5^(/,z(A)7)^(A)(7. (6.17) 

AeA 

Then [g,'y,A) is a Gabor frame, and in this case we call 7 (resp. g) the dual Gahor 
frame or simply the dual frame to g (resp. 7). 

It is shown in [22] that if S* is a Gabor frame operator, both S and the inverse frame 
operator S~^ commute with time-frequency shifts which generate S, i.e., 

ST^M^ = T^M^S and S~^T,M^ = T^M^S-^ V(x,cu)gA. (6.17) 

Thus we have 

f = SS-^f = J2{S-'f,z{X)g)z{X)g 
AeA 

= J2(f,z{X)S-'g)z{X)g (6.17) 

AeA 

where (6.17) hold since S, and hence S~^, is self-adjoint. Thus S~^g is a dual window 
for g, called the canonical dual window. We conclude that the dual frame F of a Gabor 
frame G has the same Gabor structure as G. The same property holds for S'^, i^ ^ 0, 
and in particular for S~2 , If we denote by S'^{g, A) the action of the u^^-power of the 
Gabor frame operator S on the entire Gabor frame {g, A), this property is denoted 

S''ig,A) = iS''g,A). (6.17) 

We demonstrate this in Appendix B.l. 

6.4 Lowdin orthogonalization in the continuous setting 

We describe explicitly the role that Lowdin orthogonalization plays in our OFDM 

\T 0] 
scheme. Fix a function / with good decay, let L = p t)e the generator matrix 

for {f,T,F), and assume that {f,T,F) is hnearly independent. Since {f,T,F) is a 
countable set, we can index its elements by Z; it is however more transparent to index 
the elements of (/,T, F) by Z x Z. Define {fj,k}jM^ '■= U\T,F). Since {fj,k}j,ke^ is 
linearly independent, it forms a Riesz basis for its closed linear span, and we conclude 
that there exist constants A', B' G IR+ such that 

^'||c||2< \\J2c,,kfj,k\\<B'\\ch (6.17) 

j,k£Z 
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holds for all finite sequences {cj^k}- Let J be some finite index set. Then for any finite 
sequence {cj,^}, we have 

{Rc,c)= 2_^ifj',k',fj,k)Cj',k'Cj,k=\\ 2_^ ^j',k'fj',k>\\ ■ (6-17) 

jr',fceJ j',k'eJ 

Since any cq G (?{7j) can be approximated arbitrarily closely by some finite sequence 
c, the action of R is bounded above and below on the unit ball in i (Z) by (6.4). Also 
{Re, c) = c*Rc > 0, and we conclude that (6.4) implies that the Gram matrix R is 
a positive invertible operator. From this we can infer the existence of R~2. We can 
construct an orthonormal system (0, T, F) from (/, T, F) via 



E^MAo/m, (6.17) 



k,i 

where 

Rk,e,k',u = {fk',i', fk,e)- (6-17) 

We will write 

= Lo(/,T,F) (6.17) 

for the function computed by this procedure. The proof that generates an orthonor- 
mal system is the subject of the next section. 

6.5 Lowdin orthogonalization via Gabor analysis 

We derive here the general symplectic case of Lowdin orthogonalization. 

6.5.1 Relevant theorems and lemmas of Gabor analysis 

In this section A will always be a symplectic lattice in M'^. 

Definition 6.5.1 Let /, 7 G L^(]R'^), and suppose that (/, 7,A) is a Gabor frame. 
The action of frame operator on h E L^(M'^) is 

Sf,,,^{h) = a^)J2{h, ^(A)/)^(A)7. (6.17) 

AeA 

Theorem 6.5.2 (Janssen's representation) [27] 

Sf,,,A{h) = e(A) J2 (^' ^(^°)/) ^(^°)^- (6.17) 

Aoga° 

We will also cite the following important theorems: 
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Theorem 6.5.3 (Wexler-Raz biorthogonality relations) [60] Suppose that Dg and D^ 
are bounded on i (Z*^). Then the following conditions are equivalent: 
(a) 5/,^,A = S^j,A = 1 on L\R'^) 
(b)m'{^,zi\l)f)=5x^o forKeA°. 

Corollary 6.5.4 A Gabor system (/, A) is a tight frame if and only if (/, A°) is an 
orthogonal system. In this case the frame bound A satisifes A = ^{AY- 
Proof. (From [22].) We have 

j3f,f = I (6.17) 

where / is the identity on L^(]R^). Using 7 = ^/ in Theorem 6.5.3 yields 

^-^{z{Xo)f,z{K)f)=K,K- (6.17) 

For Ao = X'o we get A = C,{AY\\f\\2 for the frame bound. Conversely, if (/, A°) is an 
orthogonal system. Theorem 6.5.2 yields Sfj = ^(A)'^||/||2/. Thus (/, A) is a tight 
frame. □ 

Theorem 6.5.5 (Ron-Shen duality principle) [45] Let f E L'^{R.'^) and a,P > 0. 
Then the Gabor system (/, A) is a frame for L^(]R°') if and only if (/, A°) is a Riesz 
basis for its closed linear span. 

It is interesting to note that this duality relation can already be extracted from a 
paper of Rieffel on non-commutative tori [42] . 

Lemma 6.5.6 Let Tio denote the closed linear span of (/, A°). // (/, A°) is a Riesz 

basis for Tio, then {Srfj^f, A°) is an orthonormal basis for Ti^. 
Proof. 

(/, A°) is a Riesz basis for TYq 
^ (/, A) is a frame for L'^{W^) (6.17) 

^ ('5/!/,a/. A) is a tight frame for L^{R'^) (6.18) 

_ 1 
^ (S'r ?^/, A°) is an orthonormal basis for Tio (6.19) 

where we used Theorem 6.5.5 for (6.17), (6.18) derives from Lemma 6.2.1, and (6.19) 
is due to Corollary 6.5.4. □ 

Lemma 6.5.7 // (/, A) is a frame for L^(]R'^), then ^(A) > 1. 

This lemma is proved in Section 9.4 of [22]. One of the consequences of Lemma 6.5.7 
and Lemma 6.5.6 is 

Corollary 6.5.8 //(/, T, F) is an orthonormal system in L^(]R'^), then TF > 1. 



54 



6.5.2 Lowdin orthogonalization and the tight frame operator 

Lemma 6.5.9 Let A be a lattice in M^"' with ^(A) > 1. Fix f E L^(]R'^), and denote 
by fx members of the Gabor system (f,A), and by fx„ members of the Gabor system 
(f, A°). Suppose (/, A°) is a Riesz basis for its closed linear span, denoted Tio- Denote 
by (i?/,A°)AoA'„ = (/a'„,/ao) ihe Gram matrix of{f,A°). Then 

^ ■■= E (^/~i)A.o /ao = S^IJ. (6.19) 

Aoga° 

Furthermore, either side of (6.19) generates an orthonormal system on A° . 

Proof. The right-hand side is in fact an orthonormal system on A°, by Lemma 6.5.6. 
It remains to prove that the sides of (6.19) are equal. For convenience we will denote 
Skjj by Si. We have, (cf. (6.14)), 



S-^z{\)f = Y.{s;--z{\)f,Sr'z{}^)f)S-^z{\')f (6.20) 

A'GA 

= Y. (<^)f ' ^r'^(^')/) s;h{x')f (6.21) 



A'eA 



since V/i G 7i, and for any frame operator S, S 2h = S°'S ^S'^h and powers of S 

_ 3 

are self-adjoint. Si * f generates a frame by (6.14). Janssen's representation of the 

_ 1 

action of the frame operator Si ^ on arbitrary h is 

s;Hh) = a^Y J2 (s;h,4Xo)s;'^f)z{x.)h. (6.21) 



AoeA° 



_d 3 



Write Ri := -R/,a°, v '■= C(A) * Si * f. If we can show that 

i?r' = R2 ■■= e(A)'i?.,A° = a^Y{{zi\o)v, ziK)v)} (6.21) 

our task will be complete. The (Ao, A'") entry of .^(A)^'^i?ii?| is 

a^r E {4>^o)f,ziK)f)J2{4K)v,ziK)v){ziK)v,ziK')v) (6.22) 
A'oeA° A'o'eA 

= e(A)^ E {zi\o)f,ziK)f){ziK)S;K,ziK')v) (6.23) 

a'oGA° 

= e(A)^(^(Ao)5i(5r^/),^(A'o» (6.24) 

= a^Y{z{x,)s;'^f,z{K')s;^f) (6.25) 

= Kx'^' (6.26) 
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where (6.23) derives from (6.21) and Proposition 2.6 in [27], (6.25) is due again to 
Proposition 2.6 in [27], and (6.26) holds by 6.5.3 on the duals in (6.14). Thus we have 
proven that 

i.e.. Lemma 6.5.9 holds by Theorem 6.5.2 with h = f. □ 

In the sequel we shall call OFDM using the above algorithm pulse-shaping OFDM, 
or PS-OFDM. 

6.6 Optimality of = Ld (/, A) 

The following theorem shows that the orthogonalization procedure described above 
yields pulse shapes that are optimal in a certain sense. 

Theorem 6.6.1 Let {fx\\ E A} with f\{t) = [z{\)f]{t) and \\f\\ = I be a Riesz 
basis for the separable Hilbert space Ti. 

(i) The function (f) = Ld (/, A) is optimally close to f in the sense that solves the 
optimization problem 

argmin \\Af — Ah\\2 when {h,A) is an orthonormal basis for Ti. 
(a) Furthermore 

\\Af - A0II < max{|l - ejf, |1 - e;±\'}, (6.25) 

where ^min (ind Omax (11"^ the minimal and maximal eigenvalues respectively of R. 

Proof. By the orthogonality relations of the ambiguity function (e.g., see Chapter 
4.2 in [22]) we have that \\Af - Ah\\ = \\f - /if. Now (6.6.1) follows immediately 
from Theorem 2.1 in [30] by using the duality relations between the systems (/, A) 
and (/,Ao). 

To show part (ii) we again use the duality relations. Recall from Section 6.5.2 that we 
can write = C,{A.)^2S~2f where S is the Gabor frame operator of the frame (/, Ao) 
(the normalization by .^(A)~2 yields ||0|| = 1). There holds ^min < C(A)S' < ^max, see 
Theorem 3.1 in [27]. The estimate (6.6.1) follows from the calculation 

ll/-^ll = ll(/-e(A)-^5-^)/||<||/-e(A)-^5-^||op 

< max{|l-l/v/^|,|l-l/v/^l}- (6.25) 
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6.7 Invariance of (p = Ld {g, T, F) under the Fourier transform 

In this section we shall treat rectangular lattices only. Recall that T is the symbol 
length and F is the carrier separation for an OFDM system. 



From here on we will always denote Gaussians by g^, i.e., 



(6.25) 



Since the pulse shape obtained from the Gaussian via (6.4) plays a major role through- 
out the paper we introduce a separate notation for this pulse to distinguish it from 
other pulses one might use for OFDM. We denote 

V?^:=Lo((?^,T,F), (6.25) 

and if cr = 1 we skip the subscripts and simply write (p or g, respectively. 





Figure 6.1: OFDM window produced by Lowdin orthogonalization (left) and its log-plot (right). 
Note the similarity to the Gaussian, except for the zero-crossings, and the exponential envelope 
decay evident in the log- plot. 



In Section 6.6 we made no mention of the TFL properties / must have for <^ to be a 
good candidate for an OFDM system. Theorem 6.6.1 says that the better the TFL of 
our Gabor atom /, the better the TFL of the Lowdin orthogonalized pulse (p. It is thus 
natural to begin with a normalized Gaussian pulse^ g. In addition to possessing the 
best possible TFL, g is invariant under the Fourier transform JF as well. This property 
is inherited by ip := Ld{g,T, F) provided that T = F and TF G Z. Indeed, for all 
/ e (L^ U L2)(M),we have {MyT^fJ = TyM_^f = e'^^'^'^M.^Tyf. Thus we see that 

^We shall show in Chapter 7 that in certain situations, g should be chirped and/or dilated. 
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under the Fourier transform, the roles of T (symbol length) and F (carrier separation) 
are rotated by ^, and a phase factor depending on the product xy is introduced. So, 
ii T = F and xy = kTiF G Z, the entire Gabor basis {g, T, F) is invariant under JF 
(since g is invariant under J-') . Or, if we have only that xy = kTiF G Z, then still ip 
is invariant under JF. In this case. 



^V = ^^Rk,lo,o9k,i (6.26) 

k,e 

= / ,Rklo o^9k/ (Fubini's Theorem) (6.27) 

k,e 

= y^ Re,k,ofl9^,k (R- is Hermitian) (6.28) 

k,e 
= ^■ 



6.8 Condition numbers and TFL 

Theorems 5.4.1 and 6.6.1 imply that the OFDM pulse shapes (p„ = Ld{g„,T,F) 
are optimally time-frequency localized among all possible OFDM pulse shapes in the 
sense that they are closest (in the L^-sense) to the Gaussian g„. Unlike g^ we see 
from Example 2 in Section 6 of [6] that ip^ has "only" exponential decay in time and 
frequency. Here the constant in the exponent depends on the condition number of 
Rg„,T,F- In Figure 7.6 we display the behavior of cond(-Rc,^,T,F) in dependence of the 
spectral efficiency {TF^^ for (TF)-^ = ^ withp = 1,2, .'. .',25. When (TF)-^ -^ 1 
we have cond{Rg^^T,F) — > cxd as a consequence of the Balian-Low Theorem and the 
functions ipa become increasingly ill-localized, cf. also Theorem 7.3 in [17]. Fortu- 
nately for typical values of {TF)~^ in OFDM such as TF G [|, |] (and a = F/T) 
the condition number of R is small. For instance for TF = 2 one can show analyt- 
ically [29] that cond(i?g^ 'J'/v^,2v^/t) = V^ and for TF = | a numerical inspection 
yields cond(i?^^ -^/v^.^v^/t) = V^- The second graph in Figure 7.6 corresponds to 
the condition number of the Gram matrix of a modified OFDM system, which we will 
introduce in Chapter 7. 

For a general theoretical investigation of the decay properties of pulse shapes gener- 
ated via (6.4) refer to [54]. Efficient numerical algorithms to compute ipo- that are 
also applicable if TF ^ N can be easily derived from the algorithms in [52]. 

It has been shown in [6] that for certain pulse shapes the so-called IOTA approach in 
[36] is equivalent to (6.4). In [36] the authors computed numerically the TFL of v? via 
(5.6) as Ar^ ■ Az/^ = 1.024- ^ (as compared to -^ for the Gaussian) and conjectured 
that this is the optimal value for all pulse shapes that create an OFDM system with 
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spectral efficiency |. Theorem 6.6.1 provides additional evidence for the correctness 
of their conjecture (as long as we consider only rectangular time-frequency lattices). 



6.9 Biorthogonal frequency-division multiplexing and OFDM- 
BFDM interpolates 

Recently biorthogonal Frequency-Division multiplexing (BFDM) has been introduced 
[32]. The reason for considering BFDM is that by giving up the orthogonality con- 
dition one can improve the TFL of the transmitter pulse shape for TF > 1. This 
comes however at the cost of increased sensitivity to AWGN, since we no longer have 
an orthogonal set of transmitter functions. Nevertheless the concept of using gen- 
eral time-frequency lattices in order to optimally adapt the transmission functions 
for time-frequency dispersive channels can be applied without any modification to 
BFDM. 

BFDM and OFDM can be seen as special cases of the following general framework 
which allows one to fine-tune the balance between AWGN sensitivity and ISI/ICI 
robustness. Assume we are given a non-orthogonal linearly independent set of trans- 
mission functions {h, A) (here h may be the Gaussian function) and let Rh,A denote 
the Gram matrix as defined in (6.3.4). 

Let R = Rg,A° ■ A pair of perfect reconstruction transmitter and receiver function sets 
(0, A) and (/, A) can be constructed by defining ([27]) 

<^ = E^MAoV, / = E^MAoV, Pe[0,l]. (6.27) 

k,e k,e 

Clearly the choice V = \ leads to a standard OFDM system (with = /) with the 
pulse shape construction as used throughout the paper. The choice p = or p = 1 
corresponds to BFDM. For < p < 1/2 (or 1/2 < p < 1 respectively) we obtain a 
transmission system that "interpolates" between OFDM and the standard setup of 
BFDM. We term these such systems extended-biorthogonal frequency- division multi- 
plexing, or EBFDM, systems. Since is the pulse whose modulations and translations 
are subject to the vagaries of the wireless channel, if has good TFL the signal is 
likely to be more robust against multipath and Doppler dispersion. The EBFDM sys- 
tems (0, A) and (/, A) become "more and more orthogonal" as p ^ | (the condition 
number of -R^^a and R/^a approaches 1), but the TFL of might degrade. Hence 
by varying a we should be able to emphasize either the AWGN performance or the 
ISI/ICI performance. 
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6.10 Numerical simulations 

For the computer simulations we assume the channel to be wide-sense stationary with 
uncorrelated scattering (WSSUS), cf. [3]. We consider a Jakes-type Doppler behavior 
and exponential decay in the temporal direction for the scattering function [41]. In 
addition we assume that tqUq < 1, i.e., the channel is underspread [3]. Following other 
work on pulse shaping OFDM [58, 24, 32, 7, 36], we consider transmission systems 
with spectral efficiency p = 0.5 (in practice a higher spectral efficiency is certainly 
preferable). 

Recall from (5.4) that the received data Ck/ can be written as 

Ck,e = ^ Cfc',/'(H0fc',r, (t>k,i) + {£, (pk,i)- (6.27) 

k',l' 

Let Su denote the time- varying transfer function (also known as Weyl symbol) of H, 
see [3] . For underspread channels and pulse shapes that are well concentrated in time 
and frequency it has been shown in [32] for rectangular time-frequency lattices that 

(H0fc/, (pk/) = S-H.{h, i^i), (Afc, i^i) e A (6.27) 

and that {ii(f)k',i', (pk/) ~ for /c 7^ /c', / 7^ /'. This result also holds for general lattices. 
We thus can use a one-tap equalizer (simple division by the value of Su at one point in 
the TF-plane) to approximately recover the data. The received data can be equalized 
by computing 

dki = ^ .''\^ , , (6.27) 

where a is the variance of the AWGN. 

Furthermore, for a fair evaluation of the robustness of CP-OFDM, PS-OFDM (and 
later LOFDM, defined in the next chapter) against ISI/ICI, we must exclude any ad- 
ditional effects resulting from coding or channel estimation. Thus we assume that the 
channel is known at the receiver and do not perform any coding/ decoding of the data. 

We measure performance by the received signal's signal-to-noise-plus-interference ra- 
tio, defined by 

SNIR(H,s) = -lOlog.o f^^) , (6.27) 

V I|C||2 / 

where c is the vector of transmitted symbols, and d is the vector of received, equalized 
(by use of (6.10)) symbols. 
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6.10.1 Our numerical implementation of OFDM and PS-OFDM 

We discuss briefly the numerical implementation of an OFDM system based on an 
orthogonalized Gaussian. First, one should have the discretized Gaussian invariant 
under the DFT, in order that the Lowdin-orthogonalized pulse ip be invariant under 
the DFT (recall Section 6.4 for the continuous case). Thus can we consistently re- 
flect the continuous signal built from ip in both the time and frequency domains, and 
avoid any issues stemming from different sampling intervals in the time and frequency 
domains. Note that if the sampling interval of the underlying continuous function 
is At, then DFT{/(nAt)} = DFT{/}(A;/(iVAt)), where N is the DFT length and 
n, k are time and frequency indices, respectively. Thus, to enforce DFT-invariance of 
discretized ip, we must have equal spacing in the time and fequency domains, and so 
we demand At = (iVAt)-^ which yields At = l/y/N. 

We choose a pulse length of A^ = 512 yielding At = 1/(16^2). Let {g, T, F) with T = 
F and g = v^e~'^* be a Gabor system with oversampling by a factor of two; i.e., T = 
F = 1/a/2, which yields a system with spectral efficiency p = | (of course in practice a 
higher spectral efficiency is desirable, but p = | was chosen for ease of programming). 
Then by the argument in Section 6.5, our OFDM system is (v?, \/2,^/2), with (p> = 
T.k,tRl}[lo9k,i with gk,i = git - kT')e'^^'^^'\ We then have T = ToAt, where Tq = 
32 = Fq are the values we use for T' and F' while programming, since the underlying 
continuous time and frequency are irrelevant on the computer. Let Nc be the number 
of subcarriers, and N^ the number of time-slots. Let c be the vector of data symbols, 
and upsample c by F to yield cp. Taking the inverse DFT of cp yields (we shall not 
complicate things by adopting or discussing the indexing or normalization conventions 
of MATLAB's FFT function here) 

cp[j] = J^ cp^^e^-^-^^o/^ = Y, c^e''^*'^* (6.27) 

e=i 1=1 

which is the sum of the symbols times the modulators since F = \f2 , t = j dt, 
and dt = 1/yN. Multiplication by the time-shifted window completes the encoding 
operation. The signal s is then subjected to a dispersion subroutine which imparts 
to s the effects of multipath delay and Doppler dispersion. Next we add noise to s, 
demodulate via the DFT, equalize using (6.10), and compute the SNIR. 

6.11 Comparison of performance of pulse-shaping OFDM 
with cycUc-prefix OFDM 

It is well-known that inserting a cyclic prefix (CP) into the OFDM symbol at each 
time-slot, while reducing the spectral efficiency (e.g., by 20 percent in IEEE 802.11a), 
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effectively precludes ISI, as long as the CP length is longer than the effective temporal 
support of the channel scattering function. CP-OFDM is of no use in combatting ICI 
resulting from frequency dispersion due to the Doppler effect from moving transmit- 
ters and/or receivers. In this section we present the results of numerical simulations 
to this effect. 

We have the freedom of several parameters in our simulations. A noise level of 10~" 
shall mean that AWGN with a variance of 10"" was added to the signal before it 
entered the receiver. A noise vector Msig of unit variance was generated in MATLAB, 
and multiples of that same Msig were used in every experiment in this section and 
Section 6.12. The simulated channel scattering function is governed by four parame- 
ters, To, z/q; 01-, and [3. The value of tq represents the maximum multipath delay (now 
in units of At), and z/q represents the maximum Doppler spread (now in units of 
A/ = At). For tq = vq < 22 the channel is underspread (the IEEE 802.11a standard 
cyclic prefix length is 20%, this is tq = t'o < 7 in our setting). The parameter a < 1 
is the maximum (attenuated) amplitude, relative to the line-of-sight strength, of the 
reflected part of the signal arriving at a delay of tq, while j3 determines the proximity 
of the maximal frequency shifts determined by vq to the asymptotes of the Jakes 
PSD (see [41]); alternatively, 13 determines how tightly clustered is the incoming sig- 
nal around zero Doppler. The channel spreading function is the scattering function 
multiplied by a normally distributed noise matrix Mscat with variance 1. Again, that 
same Mscat was used in every experiment in this section and Section 6.12. 
Note that even under the absence of the Doppler effect, (/?, optimally adapted to chan- 
nel conditions where tq = uq, usually outperforms the square pulse with cyclic prefix 
if the cyclic prefix length is smaller than tq. If however, we introduce the condition 
that z/q = To, CP-OFDM is defenseless. 

Where the Doppler effect is present, cyclic-prefix square- wave OFDM cannot mount 
any resistance to its frequency dispersion, and the resultant interference destroys the 
signal. But the signal generated from ip wards off the dispersive effects of Doppler 
interference. 



6.12 Performance of EBFDM 

In this section we consider the performance of EBFDM systems as a function of the 
exponent p. Is it well-known in communication theory that orthogonality {p = \) 
yields optimal robustness against AWGN. 

However, for nonzero dispersion the deviation from optimality of p = | for nonzero 
dispersion is often small, but thus far unexplained. Random effects are a possible 
culprit, although lengthening the signal has produced essentially the same results. 
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20 - Ricean Fading 



Ricean Fading 
cc=.50 



-*■ CP-OFDM 



(a) 



(b) 



Rayleigh Fading 



Rayleigli Fading 



(c) 



(d) 



Figure 6.2: SNIR vs. cyclic prefix length for both rapid (a — .05) and slow {a — .5) temporal decay 
of channel scattering function, under both Ricean (LOS-dominant) fading and Rayleigh fading. The 
SNIR of the Lowdin-orthogonalized Gaussian pulse (p is shown by the solid line. The noise level is 
10~^, and To = 8 but vq = 0. 



Additionally, it is not yet clear what, if any, dependence the SNIR has on a or (5, nor 
is clear the difference in character between the cusp-like maximum in Figure 6.4 and 
the apparently smooth maxima in e.g.. Figure 6.5; additionally the local concavity in 
Figure 6.4 fails to mimic the convexity in e.g.. Figure 6.5. We note that in Figures 
6.5 and 6.6, the difference in SNIR for the two noise levels above is negligible. But for 
a noise level of 10~^, the noise clearly plays a significant role in lowering the values 
of the SNIR. 
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Figure 6.3: The same parameters as in Figure 6.2, except h'o = tq = 8. 
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0.2 0.3 0.4 0.5 0.6 0.7 



Figure 6.4: SNIR vs. p for noise level 10 ^® and tq = uq = 0. 
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Figure 6.5: SNIR vs. p for noise level 10~^ and tq = t'o = 4 (Figures 6.5 (a) and (b)) and 
To — 1^0 — 8 (Figures 6.5 (c) and (d)) for both rapid {a = .05) and slow (a = .5) temporal decay 
of channel scattering function, and wide (Figures 6.5 (a) and (c)) and narrow (Figures 6.5 (b) and 
(d)) distributions of Doppler shifts. 
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Figure 6.6: The same parameters as for Figure 6.5 except that the noise level is 10 
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Figure 6.7: The same parameters as for Figures 6.6 and 6.5 except that the noise level is 10 
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Chapter 7 

A General Framework for Pulse 
Shape Design 



7.1 Adaptation of the pulse to channel conditions 

It has often been suggested that the time duration and frequency duration of the 
transmitter functions should be adapted to the channel properties, described by the 
channel scattering function Cjj{x, u). The ambiguity function Aipj^^ilx, u) is the inner 
product of a receiver function centered at {kT, iF) and a channel-scattered transmit- 
ter function centered at {kT±x, iF±uj). Thus Aip^^i describes the ability of an OFDM 
system to detect data symbols in the presence of time- and frequency-dispersion. For 
a function / and e > define the effective support of the ambiguity function Af by 



suppMf) = {it,uj):\Af{t,cu)\>e}. 



(7.0) 




Figure 7.1: Ambiguity function poorly matched to channel scattering function. 



Consider Figure 7.1. In region Y, {Aipk/l > e, but the channel scatters no input to Y 
so designing a receiver that can detect a F-scattered pulse is wasteful. Conversely, if 
the receiver cannot detect a pulse scattered to the region X, data will be lost. Thus 
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we should adapt the pulse so that supp^Aip^^g coincides, at least roughly, to the sup- 
port of the channel scattering function. 

In order to reduce ISI and ICI it often has been suggested to choose the time du- 
ration Atj and the frequency duration Az/j of the transmission pulse as well as the 
parameters F and T such that they are adapted to the channel properties, cf. e.g. 
[36, 24, 32] and for a more theoretical approach see [31]. This means that for given 
maximum multipath delay tq and maximum Doppler spread uq we choose Arj, Az/j, 
T, and F such that 

z/q Auf r 

For the Gaussian-based window ipcr this means that for given spectral efficiency p, we 
set o" ~ ^ and T = \/ -^, F = ./-. Using basic properties of the dilation operator 
Va we obtain the relation 

Va = Ld {g^, 1/v^, \/o7p) = V^{Ld {g, 1/^, l/y/p)) = V^^p. (7.0) 

A detailed proof of (7.1) is in Appendix B.2. 

Thus the transmission pulse ip^- that is matched to the channel parameters Tq and 
pq is obtained from the "standard" pulse ip by an appropriate dilation. Hence at 
least theoretically it is very easy to adapt the OFDM system to time-varying channel 
conditions. We only have to know how the ratio of Tq and Uq changes and adapt the 
lattice and (p^- accordingly by a simple dilation. Of course, in practice one must also 
deal with issues such as time and frequency synchronization when one changes the 
lattice parameters, which may be a non-trivial task. 

7.2 Lattice-OFDM, wireless channels, and interference 

The robustness of an OFDM system against ISI and ICI depends essentially on two 
factors: 

(i) the TFL of the functions ipki] 

(ii) the distance between adjacent lattice points {kT,iF), {{k ± 1)T, ((/ ± 1)-F). 

We have seen in Section 7.1 how to construct for given channel conditions OFDM 
pulse shapes with optimal time-frequency concentration. If we want to further im- 
prove the stability of the OFDM system against interference we have to increase the 
distance between adjacent data symbols by increasing the distance between adjacent 
grid points. However increasing T and/or F results in an undesirable loss of spectral 
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efficiency. Can we increase tlie distance between adjacent grid points witliout reduc- 
ing tlie spectral efficiency? 

We will show in this section that this is indeed possible. The solution lies in the 
introduction of a more general and flexible OFDM setup which we call Lattice- OFDM 
(LOFDM). 

7.2.1 Lattice-OFDM 

The generator matrix L^ for the rectangular lattice A^ displayed in Figure 7.2 is 



^R 



T 
F 



(7.0) 



the density of which is ^{Ar) = l/det(L«) = 1/(TF). 



The hexagonal lattice Ah (shown in Figure 7.3 with T = F) with generator matrix 



^H 



V2rp y/2 rp' 

fF_ 



(7.0) 



has the same density, viz. ^{Ah) = l/det(L//) = 1/{TF). 



Definition 7.2.1 Let A be a lattice with generator matrix 



X y 
z 



x,z j^ 0, 



(7.0) 



and let E L'^(R.) be given. As in (6. 3. 3), we denote by (0, A) the function system 

0H(t) = 0(t - Afe)e=^"'^^*, (Afe,/i,)GA, keZ,i = 0,...,N-loTieZ. (7.0) 

If the functions (p^e are mutually orthonormal we call (0, A) a Lattice-OFDM (LOFDM) 
system. ""^ 

It is clear that the spectral efficiency p of an LOFDM system (0, A) coincides with 
the density of the lattice; i.e., p = l/det(L). 



^We have restricted in the definition of LOFDM to lattices with upper-triangular generator matrices, instead of 
general invertible 2x2 generator matrices. The reason is that unlike general lattices those with upper-triangular 
generator matrix - though not necessarily separable - still possess a distinctive orientation along the time-axis and the 
frequency- axis respectively. This orientation is very useful for an efficient numerical implementation (see also Section 
7.3) and it keeps the time delay introduced by the non-separability of the lattice small. All results that follow could 
be derived (with some minor adaptations) for general lattices, however with little or no practical benefit. 
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7.2.2 Sphere-packings and minimization of ISI/ICI 

In this section we show how to design pulse shapes and lattices which provide opti- 
mal robustness for time-frequency dispersive channels. We start with a specific but 
important example. The analysis of the general case will be based on this example. 

Assume that the parameters of the channel scattering function satisfy tq = uq (we 
will discuss more general choices for ro,z/o later). By the results in Section 6.6 the 
optimal pulse shape ip in this case is given by v^ = Ld{g,T,F) with g{t) = 2'4,e~'^* 
(since a = ^ = 1) and T = F. 

The ambiguity function oi g is rotation-invariant, see [18], therefore 1^(7^/1 is rotation 
invariant. Hence for any e > 0, supp^{Agk/) is a sphere with center {kT,iF). It 
follows from (6.6.1) that if the matrix R is well-conditioned then Aip ^ Ag. Hence 
supp^{Aipk,e) is approximately a sphere with center {kT,£F). The OFDM system 
{(p, T, F) with T = F and a = 1 can thus be represented schematically by Figure 7.2. 
The data symbol Ck,e is transmitted at the lattice point {kT,iF). The gray spheres 
represent supp^Aipk/- We have chosen e such that the difference between the setup 
in Figure 7.2 and Figure 7.3, see below, becomes obvious. Certainly the spheres may 
overlap for smaller e, but the A(pk,e are mutually orthogonal, since the (pk/ are mu- 
tually orthogonal. However after passing through the channel H the orthogonality 
between the ipk,e (and thus between the Aipk/) is lost and the support of Aipk/ is 
spread out in time and frequency. And now the overlap between adjacent Aipk/^ 
(or rather between AHipk/ and AHipk'^u) plays a crucial role. Roughly speaking less 
overlap between the spheres results in less ISI and ICI. 

Note that since 'Pk/{t) = 0(e~'*') we have 

I {iiipk,e, Vk±i,i±i) I > I (Hv9fc,^, Vk±n,i±n) I , forn > 1. (7.0) 

Therefore we can concentrate on the interference between Ck/ and direct neighbor 
symbols Cfc-ti^^-i-i, since the interference of Ck/ with symbols Ck±n,i±n where n > 1 will 
be negligibly small by comparison. Hence the problem reduces to maximizing the 
distance between immediate neighbor lattice points (A;, £) and (/c ± 1,/ ± 1) without 
losing spectral efficiency. The answer to this problem comes from the theory of sphere 
packings. 

The sphere packing problem can be loosely described as follows. We are given an 
infinite number of rf-dimensional spheres of equal size. We identify the midpoint of 
each sphere with a point of a rf-dimensional lattice A in W^. In the sphere packing 
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(a) 



(b) 



Figure 7.2: An OFDM system with rectangular lattice in the time-frequency plane. The data 
symbols Ck,i are transmitted at the lattice points {kT,£F). The spheres represent the effective 
support supp^Aipkj of the ambiguity function Aipk^i- The choice of e is such that the difference to 
the LOFDM setup displayed in Figure 7.3 becomes obvious. 



problem one tries to find the lattice Aopt that solves 

f Volume of a sphere 



max 

A 



V 



Density of A J 

It is easy to see that in the sphere packing problem the radius r of a sphere is given 
by 

{l|Afc-A,||}). (7.0) 



mm 

2 VAfe.Ai e A, 

In words, the radius is equal to half of the minimal distance between two lattice 
points. We can reformulate (7.2.2) as follows: 



max( min {||Afc-Az||}) 

A Afc,A,eA, 



subject to ^(A) = ^0 



(7.0) 



for some arbitrary ^o > 0. Obviously we find ourselves confronted with the same 
problem as in OFDM, when we want to maximize the distance between adjacent lat- 
tice points in order to minimize ISI/ICI while keeping the spectral efficiency constant. 



It is well-known that the sphere packing problem in two dimensions is solved by the 
hexagonal lattice, cf. [12]. Thus for the case tq = uq an OFDM system constructed 
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from a hexagonal lattice as shown in Figure 7.3 will provide optimal protection against 
ISI/ICI due to the increased distance between adjacent lattice points. The effect of 
the increased distance becomes apparent when comparing Figure 7.2 with Figure 7.3. 
In both cases we have chosen the same e as the threshold for the effective support of 
the ambiguity function. 



1 /""^ 


"o 


= ^0 






/ 
/ 




T 




(a) 



(b) 



Figure 7.3: Optimal LOFDM system (Fig. 7.3(b)) for the channel scattering shown in Fig. 7.3(a). 
The data symbols Ck,i are transmitted at hexagonal lattice points. The spheres represent supp^Aipk^i 
with the same e as in Figure 7.2. The use of the hexagonal lattice clearly increases the distance 
between adjacent spheres as compared to the rectangular lattice (cf. Figure 7.2), thus leading to 
reduced ISI/ICI for time- frequency dispersive channels. 



Furthermore when we construct an orthogonal set of transmission functions (y?, A) 
from a non-orthogonal set of Gaussians {ga,h) via ^o- = Ld{gcr,A) the increased 
distance between adjacent lattice points yields a smaller correlation |(v^o-a'' V^o-a)I be- 
tween neighbor functions than in the case of a rectangular lattice. This results in 
a smaller condition number of the associated Gram matrix Rg^^A for Rg^^T,F at the 
same spectral efficiency, as can be seen in Figure 7.6. 

Of course the assumption tq = vq in our example is very restrictive. In the next 
subsection we consider more general channel conditions and derive a general LOFDM 
design rule in order to minimize the joint ISI/ICI. 



We shall from this point forward restrict ourselves to d = 1 since we are constrained 
to use only symplectic lattices. In M? any lattice is symplectic, but this is not the 
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case in M with d > 2. 

7.2.3 Optimal LOFDM design for the minimization of ISI/ICI 
Definition 7.2.2 We define the dilation operator "Dq, by 

VJXt) = V^fiat), a e M+, (7.0) 

and multiplication by a chirp is 

Cpf{t) = f{t)e--^^'\ /3gM. (7.0) 

Va and Ci3 are unitary linear operators on i7^(M), and we have V'^^ = V^-i and 



^/3 - C-/3- 




For f,heL'^{R) 


we have [22 




A{VJ,V^h){t,uj) = A{f,h){at,a-'uj) 




AiCpf,Cph)it,iu) = Aif,h)it,u; + (3t). 


Lemma 7.2.3 






J^V^-i = V^J" 


Proof. 










= ^ / /(t/a)e-2™*rft 








= V^J^fit). 



(7.0) 



Recall that tq and z/q describe the support of the channel scattering function. More 
precisely, assuming that the support of the scattering function is circumscribed by a 
rectangle (or ellipse), then tq gives the length of the rectangle (in the time-direction) 
and z/q corresponds to the width (in the frequency direction). In the following the- 
orem we extend the result from Section 7.2.2 to more general (and more realistic) 
choices for tq and vq. The theorem is essentially only a special case of a deep the- 
orem by Grochenig [22], but it shows how we can switch between different lattices 
by adapting the pulse shape accordingly in order to preserve the properties of the 
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optimal hexagonal case. In fact our approach will also include the case where the 
support of the scattering function may be concentrated along a rotated rectangle in 
the time-frequency domain (like in the case of a chirp). Such a setup has also been 
investigated in [1]. 



Theorem 7.2.4 Let (/i,Ai) be a non-orthogonal linearly independent system with 
||-R/i,Aill < oo, and denote ipi = Lo(/i,Ai). Assume A2 is a lattice with density 
^(A^) = ^(A^) and let 



Li 



xi yi 

Zl 



and Ln 



X2 1/2 
Z2 



he the generator matrices of Ai and A2, respectively. If 
with 



^2 n y2 . VlZi 

" = — , P = + 2' 

Zl Z2 Zl Z2 



then 



and if we denote ip2 := Ld (/2, A2), then 

(P2 = J^Va-iC_i3a'2j^~^(Pl 



(7.-4) 

(7.-4) 
(7.-4) 
(7.-4) 

(7.-4) 



Proof. The proof consists of applying Proposition 9.4.4 in [22] twice. Fix an arbitrary 
rectangular lattice Af> with generator matrix 



'R 



T 
F 



(7.-4) 



such that ^(A/j) = 'C(^i)- It is known that any A G Sp((i) can be written as a finite 
product of matrices of the form J, Cp, and D^- It is thus clear with some calculation 
that Li and L2 can be expressed as 



J 



■^1 = 


JCp^Da^J 


-'Lr, 


with «! = 


Zl 


Zl 


(7.-3) 


J2 = J C 13^2^02 J Lr, 


with 02 = 


_ ^2 n _ y2 

-r^'-~V2 


(7.-2) 


are defined by 








' 1" 
-1 


! Cl3 = 


"1 0" 

p 1 


, D^ = 


a 



" 


, a,/3 G M. 


(7.-2) 
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Let Sp(l,]R) denote the symplectic group^, and let /i(-) denote the symplectic (uni- 
tary) operator appearing in the metaplectic representation of Sp(l,]R) (see (B.3.1) in 
Appendix B.3). Let J^Cp^Da be defined as in (7.2.3), then 

^(j)=^-i ^{Cp)=C^p i2{D^)=V^-i (7.-2) 

Furthermore there holds (despite the fact that fi is not a homomorphism in general) 

(7.-2) 



Define 
Then 



l^iJCfsD^J-') = /i(J)^(C^)^(D„)/i(J-i). 
Ai := JCfs^Da^J , A2 := JCp^Da^J ■ 

= J C 13^0^2 -13^ J 



(7.-2) 



JDc,2C 



J-\ 



(7.-4) 



where in the last step we used the fact that CfiDa = DaC/sa^. From (7.-4), we conclude 



/x(A2A^^) = fi{JD^C 



ft(^) -ft 



J-'] 



J^~^ V^ c 



ft-ft(^) 



ft(^) -ft 

21 J^ 



)KJ 



-i\ 



= J^~^ V^-i C_i3c,2 J^. (7.-6) 

Set fji = fi{Ai)^^f\ and /2 = ^4(^2) ""^/ij- From Proposition 9.4.4 in [22] we have that 

(/H,A^)=/i(^i)"'(/i,Ai), (7.-6) 

and 



(/2,A2)=M^2)(/iJ,Aij). 

Since /i(Ai) and fJ.{A2) are unitary, we get 



-R/i,Ai - Rfa, 



R 



hM 



(7.-6) 



(7.-6) 



where L2 = A2A~['^Li. Equation (7.2.4) now follows from (7.-6), (7.2.4) and (7.2.4). 



^Any lattice in R^ has a generator matrix in Sp{l, 
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It is possible to prove Theorem 7.2.4 directly using signal analytic methods. However 
we feel that the metaplectic representation is simpler and more elegant. The reader 
familiar with the mathematical theory of quantum physics will notice that Theorem 
7.2.4 is essentially a consequence of the celebrated Stone-von Neumann Theorem on 
unitary representations of the Heisenberg group (see Appendix B.3). 

The main feature of Theorem 7.2.4 is that it yields a simple rule for the design of 
optimal LOFDM systems for time-frequency dispersive channels. We illustrate this 
by several examples. 

7.2.4 Some examples 

Example 7.2.5 We compute (/9/j, the pulse generating a tight frame based on a hexag- 
onal lattice with twofold oversampling, from (^^, the pulse generating a tight frame 
based on a rectangular lattice with twofold oversampling (these pulses could then be 
used as generators of OFDM systems on their respective adjoint lattices). If 



2 ^ 

^ 

^ 2 



and Lh = ^ 



■1 \ 

4, 



then a = f /f = f 
to get 



and P 



-y=, and we simply use (7.2.3), (7.2.3), and (7.2.4) 



^h 



J^V^ClT ^(Pr- 



(7.-7) 



n 



In Section 6.6 we mentioned that the TFL of the pulse ip = Ld {g, \/2, \/2) measured 
in terms of second-order moments is At^-Au^ = 1.024~, which is conjectured in [36] 
to be the optimum among all pulses creating an OFDM system (or OQAM-OFDM 
system) . This conjecture is plausible for rectangular lattices, but is untrue if we allow 
more general lattices. For instance if we use the hexagonal lattice A^ as defined in 
(7.2.1) with T = F = v2 and set ipn = Ld {g, Ah) we compute numerically the TFL 
Ar„„ ■ Au, 



'fH 



H>H 



1-014 -4^. 



Conjecture 7.2.6 ipu yields the optimal TFL among all possible lattices and all 
possible pulse shapes that create an LOFDM system, with spectral efficiency p = 1/2 
(and similarly for other values of p). 

In the following examples we consider LOFDM systems with spectral efficiency p. For 
the remainder of this section Ai will denote the hexagonal lattice A^ with generator 
matrix 






n 





2v/3 
V2 . 



(7.-7) 
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and 



ifi := Ld{g,Ah). 



(7.-7) 



We have seen in Subsection 7.2.2 that the LOFDM system ((/9i, A^) is optimally ro- 
bust against ISI/ICI if tq = vq. Based on Theorem 7.2.4 we are able now to consider 
more general situations. 



Example 7.2.7 Fix p and assume Tq = -z/q for some ct > 0. The LOFDM system 
that is optimally adapted to these channel conditions with respect to minimizing 
the joint ISI/ICI is given by {(pa-,-Aa)) where the generator matrix L^j of Ag- can be 
computed by left-multiplication of Li by the dilation matrix -Di/^ as ^ 






T T/2 
F 



with T 



V2 



'=T 



(7.7) 



and the optimal pulse shape is ip„ = Ld {g„, K„). 

Consider for instance tq = \/3z/o and p = \. By using (7.2.7) we obtain 



2 1 
1 



(7.-7) 



as generator matrix of the optimal lattice A2. Using (7.2.4) and (7.2.4) with Li as in 
(7.2.4) and ifi as in (7.2.4) we compute a = i, P = 0, hence 






(7.-7) 



where we have used Lemma 7.2.3. Thus the optimal LOFDM setup for the given 
conditions p = |, Tq = VSuq can be schematically represented as shown in Figure 7.4. 

Example 7.2.8 Assume p = | and tq = vSt'o? but now let the scattering function 
be rotated clockwise by 45 degrees in the time-frequency plane, as shown in Figure 
7.5(a). In this case Theorem 7.2.4 together with the adaptation rule (7.1) yields 



L'2 = Rotri-Ab) ■ Li 



/2 V2] 




[2 ll 




2 2 









V2 V2 




n 1 




2 2 J 









V2 V2 
-V2 



(7.-7) 



where i?otr(— 45) is the left-multiplication Givens rotator. Now L2 generates the same 
lattice as 



U-L' 



-1 

1 1 



V2 V2] 




'V2 " 


-V2 _ 




. A 



(7.-7) 



^see (7.1) and (7.2.3); also note that since it is the column vectors oi the generator matrices which are the generating 
vectors, left-multiplication is the correct operation 
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12 3 4 5 

Time 



(a) 



(b) 



Figure 7.4: Optimal LOFDM setup (Fig. 7.4(b)) for a wireless channel with tq = VSi^o (scattering 
function shown in Fig. 7.4(a)). The ellipses represent supp^A{Lpr^)k/, where we have chosen the 
same e as in Figure 7.3. 



since U is unitary and leaves one lattice-generating vector unchanged, and re-forms 
the second as a linear combination of the lattice-generating vectors in Lg. Thus we 
can take L2 as the generator matrix of the optimal lattice A2, which is a rectangular 
lattice in this case. Proceeding in the same way as in the previous example we 
calculate a = -^, /3 



|, and 



'^'^it)= TVa~iC_p^2T ViW 



V2 ^3 



(7.-7) 



Since if>r- is the Fourier transform of a dilated chirped Gaussian, formulas (7.2.3), 
(7.2.3), and (6.6.1) then show that su'p'p^Ai^)' r- is (approximately) an ellipse, with 

ratio of semi-major vs semi-minor axis equal to y^, rotated clockwise by 45 degrees. 
Hence the corresponding LOFDM system (v^'/o , A/j) can be schematically represented 
as in Figure 7.5. 

Remark: The pulse-shaping framework derived in Section 6.6 also holds for general 
lattices, since Theorem 6.6.1 carries over to general lattices (see Appendix B.3.3). 
The proof of the extension of Theorem 6.6.1 can be done by using Theorem 7.2.4 and 
the results derived in Chapter 9.4 in [22]. 



It is interesting to compare the condition number of R^ := Rg^k^ with that of Rh := 
Rg,ku for T = F = v^, where the generator matrices for these lattices are as defined 



in (7.2.1) and (7.2.1). We get cond(i?H) 
suspiciously close to \pl. 



a/2 and cond(i?//) ^ 1.2599 which is 



Conjecture 7.2.9 R, 



■3,A 



\pl is the minimal condition number achievable for any 
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^n=3-v„ 




(a) 



(b) 



Figure 7.5: Optimal LOFDM setup (Fig. 7.5(b)) for a wireless channel with 45-degree rotated 
scattering function satisfying tq = vS^'o (see Fig. 7.5(a)). The ellipses represent supp^A{i/j2)k,e, 
where we have chosen the same e as in Figure 7.3. 



lattice A with ^(A) = 1/2. 

If j/q ~ then OFDM systems employing (almost) rectangular pulses such as cyclic 
prefix OFDM are of course optimal. The essential support of the ambiguity function 
of a rectangular pulse has rectangular shape. It is clear that the optimal lattice for 
such pulses is a rectangular lattice. Thus the introduction of general lattices would 
only decrease the performance in such a situation. 



7.3 Efficiency considerations in the numerical implementa- 
tion of LOFDM 

Before we proceed to a performance comparison between OFDM and LOFDM we 
briefly discuss the numerical implementation of LOFDM. We do not go into detail 
here, but note that the results in this section can be readily derived by using the link 
between OFDM and Gabor analysis and the efficient algorithms that have been de- 
veloped in association with suitable finite-dimensional models for Gabor analysis [52]. 



The requirement for an efficient implementation places some limits on the possible 
values of the entries x, y, z of the generator matrix of the lattice (see Equation 7.2.1). 
One can easily see that an FFT-based implementation is only possible if x,y,z & Q. 
In particular if x/y G N then the lattice can be written as sum of x/y rectangular 
lattices, which differ from each other only by the time-frequency offset. For each of 
these rectangular sub-lattices we can use the usual FFT-based implementation for 
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pulse-shaping OFDM systems. An example of a "good" generator matrix is 

L 



T T/2 
F 



with TF e Q. Specifically, if TF = p/q with p, 5' € N, the length of the FFT is given 
by s'^pq, where s is a constant that depends on the ratio of the sampling rate and the 
decimation factor in the A/D conversion and is often of the form 5 = 2". 

7.4 Comparison of performance of LOFDM vs OFDM 

We compare the proposed LOFDM system to a standard pulse-shaping OFDM sys- 
tem of (5.1) associated with the rectangular lattice. For both systems we use the 
pulse obtained from the Gaussian via Ld {g„, Kr) and Lb (g^^, Ah) respectively. For 
the rectangular lattice the OFDM system obtained in that way coincides with the 
setup suggested in [36]. In both cases we optimally adapt the pulse and the lattice 
parameters to the given delay spread and Doppler spread. 

All parameters are as described in Section 6.10.1, except that AWGN on the order 
of 0.2% of the signal energy is added during the simulated transmission (the noise 
level is 9). We choose this small percentage of noise in our comparison in order to 
concentrate on the effects caused by ISI and ICI. 

For both OFDM and LOFDM we compute the SNIR of the received data dk/, which 
have been equalized according to (6.10). Since both OFDM and LOFDM are orthog- 
onal systems they are equally well (in fact optimally) adapted to AWGN. Hence any 
statistical difference in the SNIR can be attributed to the difference in robustness 
against ISI/ICI. 

We fix the product Tq ■ z/q, randomly vary the parameters of the scattering for each 
experiment and compute the average SNIR for OFDM and LOFDM over 100 such 
experiments. We repeat this simulation for a variety of values Tq • Uq in the range 
of 10~^ — 10~^ (but we always adapt both systems approximately optimal to the 
new values of Tq and z/q). The result is shown in Figure 7.7. As expected from the 
theoretical analysis LOFDM yields indeed a better SNIR than OFDM. On the average 
the improvement of LOFDM compared to OFDM is in the range of IdB. 
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Spectral efficiency 

Figure 7.6: Comparison of the condition number of the Gram matrix of the OFDM-type function 
system (g^ , T, F) (g^ is the Gaussian) for the rectangular lattice and the hexagonal lattice when 
the product {TF)~^ approaches 1. In both cases the condition number is small for TF £ [i, |]. As 
expected cond(i?) for the hexagonal lattice is definitely smaller than cond(i?) for the rectangular 
lattice, thus demonstrating the potential advantage of designing OFDM systems for general time- 
frequency lattices. 



Comparison of LOFDIUI and OFDM 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 

varying decay parameter for channel scattering function 



Figure 7.7: Comparison of SNIR for OFDM and LOFDM for tq ■ i^o ^ [10"'*, 10"^] 
LOFDM scheme yields an improvement of about IdB. 



The proposed 
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Chapter 8 

Future Research 



8.1 Weighted (under matrix and Hadamard multiplication) 
Lowdin orthogonalization via the SVD 

Let $ G C'"^" with m > n, with rank(<l>) = n. We consider the general case of 
possibly nonsquare, rank-deficient weight matrix W. Let W G C"^'^, k > n with 
rank(W^) = r < n. The (reduced) singular value decomposition, or SVD, of ^W is 



^W = UDV* 

Dl Okx{k-r) 



[Ui,U^. 



2 



U(fc— r)xr 'J{k—r)x(k—r) 



[Vl,V2[ 



= U,D,V* (8.-1) 

where U G C™'^'^, V G C^^^, D G M'^^^; Di is the diagonal r x r matrix containing 
the r nonzero singular values of $1V, f/i G C™'^'', and Vi G C''^'". The SVD of W is 
W = UwDwV^ = UwiDwiV^i, where Uw e C"^^DvK e M*^^^ Viy G C*^^^ and 
?7vKi G C"^'', DvFi G M"^", and Vwi G C^^^ 

Lemma 8.1.1 Among all unitary matrices Q, ||($ — Q)W^||ir is minimized by 

Qopt = P:= UV\ (8.-1) 

where U and V are unitary matrices in USV*, the reduced, compressed SVD of 
^WW* (U G C"^% S G M''^^ and V G C"^''). 

Proof. We have 

||<l>W^-gW^||^ = {^W -QW,^W -QW) 

= \\<^W\\l + ligW^II^ - {^W, QW) - {QW, ^W) 

= \\(!>W\\l+\\W\\l-2Re{^W,QW). (8.-2) 
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Denote the singular values in S by o"j, and V*Q*U by {tjj}. Clearly we want to 
maximize Re{^W,QW). But 



Re{<^W,QW) = Re{tT{<^WW*Q*)} 
= Re{tT{USV*Q*)} 



Re{tr{SV*Q*U) 

Re \ ^a,t,, I . (8.-4) 



i=l 



This sum is maximized when tu = 1 for all i, or when V*Q*U = Ir, which yields 
Q = UV*. □ 

The goal of the above treatment is to enable us to say something about Lowdin or- 
thogonalization in the case of Gabor frames under a weighted norm. That is, in the 
Banach space L^, where s is the weight function, which orthogonalization of the Ga- 
bor frame {g, a, j3) retains a maximal proportion of the excellent TF-localization of gl 



Obviously, the result above shows that Lowdin orthogonalization is optimal in the Ba- 
nach space described by the norm || ■ Vr||F where W is any weight factor in the matrix 
multiplication. We now formulate the discretized problem. Write gnm '■= MnaTmpg, 
the continuous translated modulated Gaussian, and write ^„„ := MnaTm/sg for the 
discretized frame element, say a row vector for concreteness. Our linearly inde- 
pendent set $ is the matrix whose rows are the gmnS. Note that (for the cur- 
rent treatment, at least) the weight function is constant for a given m, even as n 
varies; but as m varies, the weight function must translate along with the gaus- 
sian. Thus different translations of the weight function will multiply, point-by point, 
different row-blocks of $. This is a highly nonlinear operation, and thus the above 
treatment does not apply in its current form. We shall denote point-by-point ma- 
trix multiplication by and the inverse operation, point-by-point matrix division, 
by 0. Note that both and conform to the (right or left) distributive rule: 
{A + B)qC = AqC + BqC,{A + B)0C = A0C + B0C 

Let $ be as above, but the form of W is now very specific. Suppose we have a weight 
function w : M. ^ M. : t h^ w(t) with ■U7(t) > 0, Vt G M; a typical weight function is in 
the standard class of polynomial type Ws{t) = (1 + 1^)'*. Denote the discretization of 
w by w°[(i]. Then the weight matrix W will be of the form 
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w° 



w°[n-iiT] 



w°[d-i2T] 



1.-5) 



W°[d-iMT] 

where T is the symbol length, M is the total number of translations in the disrete 
model, and each of the rows w°[d — ijT] is a particular translation of w°[d] and is 
repeated N times, where A^ is the finite number of modulations in the model. Note 
that W is quite rank- deficient, but more importantly, has no zero entries (in fact for 
the polynomial weights we have Wij > 1 V i, j). 



Let $ be as above. The goal is to find, among all unitary matrices Q, the matrix Qmin 
which minimizes || ($ — Q) (DW\\f- The method used in the proof of Lemma 8. LI fails 
in the Hadamard multiplication case. I will attempt to find a method which yields 
the optimal Q for this case. 

8.2 Is C^ a syminetric Banach algebra? 

In Section 3.3.2, we noted that it seems as yet unknown whether or not C^ is a 
symmetric Banach algebra. 

Conjecture 8.2.1 Let h E L^iW'-). Then there exists ts{x) = (1 + |a;|)'' , 6 > such 
thathL G Cl^. 

If Conjecture 8.2.1 can be proven, we might conclude that C^ is symmetric, provided 
that we can redefine C^ as the union |J £^ and assign a norm to each element of 

Tg , 5>0 

h E C^ as \\h\\ci = max \\h\\c-i ■ In this case the hypotheses of Barnes' Lemma might 

5 ^^ 

be satisfied for C^. 



8.3 Extending the results regarding Banach algebras to quo- 
tient spaces 

It is possible to extend the results in Chapter 3 to Banach algebras of distributions 
containing non-invertible operators, by forming equivalence classes. Fornasier and 
Grochenig demonstrate this concept for pseudo-inverses in [19], but it is, I think, 
possible to extend their result to the more general setting of quotient spaces, based 
on the ideas of Barnes in [2]. Let £^o and £„o denote the extensions of Cl and £„ 
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containing also non-invertible operators, and let -ft'^io and K^^o denote the kernels of 
£jo and £„o respectively. 

Conjecture 8.3.1 The spaces C\ = C\ o / K^o and C'^ = C^ o / K^^o possess the 
same decay properties as Cl and C^, respectively. 

8.4 Numerical solution of integral equations and the finite 
section method 

It is often of interest to solve Ax = b where A is bi-infinite and x, b are infinite. Of 
course closed-form solutions are usually impossible, so we turn to the computer for 
our approximate solution. For any sequence y, let 

PnV = {■■-, 0, 0, y_n, y-n+l, • • • , Vn-l, Vn, 0, 0, . . .). (8.-5) 

A typical approach to finding an approximate solution to Ax = b, termed the finite 
section method (see e.g., [55]), is to truncate A via the projection PnAPn, which yields 
a finite section of A padded with zeros outside the section defined by n, and similarly 
with Pnb. We then solve the system PnAPnX = Pnb, where in general x ^ PnX. 

The rate of convergence of £ ^ a; is of interest, whether in the discrete setting 
of matrix equations or the continuous setting of integral equations. The results of 
Chapter 3 will have applicability to this question, since the kernels in Cl and C^ are 
defined by weights which by construction allow us to quantify convergence rates of 
numerical methods. 

8.5 Dependence of EBFDM SNIR on a , (3 and p 

The dependence of EBFDM SNIR on a , f3 and p is not obvious (at least to the 
author), and warrants further investigation, as does the transition from the cusp- 
like maximum at p = 1/2 for no dispersion to the smooth maximum for dispersive 
channels. 
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Appendix A 

Integral Operators and Kernel 
Decay 

A.l Associativity of • 

Let f -k g = d, g -k h = p, {f * g) -k h = q, f -k (g -k h) = r. The task is to show that 
q = r. We have 



d{x,z)= I f{t,z)g{x,t)dt and p{y,Lu)= / g{t,uj) h{y,t) dt 

Thus 

q{y,z) = I d{x, z) h{y,x) dx 

f{t,z)g{x,t)dt\ h{y,x)dx 
f{t, z)g{x, t) h{y, x) dt dx 
f{t,^)[ / g{x,t)h{y,x)dt) dx 



and thus 



f{t,z)p{y,t)dx 
r{y,z) 

{fkg)^h = fk{g^h). (A.-7) 



A. 2 Schur's test 

Lemma A. 2.1 Suppose that f is a measurable function on M^*^ which satifies 

sup / \f{t,y)\dt<Ci, and 

sup f\fit,y)\dy<C2, (A.-7) 

Then the integral operator f defined by fh{y) = J^f{t,y)h{t) dt is bounded from 
L (M*^) to L (M.'^) for 1 < p < oo, and the operator norm of f is bounded by 

||/||lp(rHlp(r) < Ci C2 (A.-7) 

where p' is the Holder conjugate of p. 
Proof. Since 1/p' + 1/p = 1, we can apply Holder's inequality to 

/ f''^'\t,y)f''%t,y)h{t)dt 

jR'i 

to get 

\fh{y)\ < (I \f'/^'{t,y)fdt)' i I \f'/^it,yr\h{t)\^dt 

< ( [ \f{t,y)\dtY( [\f{t,y)\\h{t)\^dt 
which yields 

WfhWl < I [ I \fit,y)\dty [ I \f{t,y)\\hitWdt]dy 

Jr VJRd / \Jr 

< Cf [ [ \f{t,y)\\h{t)rdtdy 

jRd- jRd- 



< Cr / / \fit,y)\\hit)\^dydt (Fubini) 

jRd- jRd- 

< Cf [ ( [ \fit,y)\dy)\hit)rdt 

jRd- XJRd J 

< C'l^'' C2\\h\\l (A.-12) 
and by taking p*^ roots, we are done. □ 
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A. 3 Properties of admissible weights 

A. 3.1 Submultiplicativity of concave k with ^(0) = 

Recall that k denotes the exponent in v{x) = e'*'-'^'-'. Since /t(0) = 0, by a simple 
corollary to Lemma 5.16 in [46] we have -j^ < -^^ =^ ^^y^ < ^(a), and 

/t(a + b) nip) 
a + b ~ b 

b 

= K{b) + %{b) <K{b) + K{a), (A.-13) 

and thus (3.3) holds. 

The proof that concave k, with k, > and k,{0) = is nondecreasing is menial 
but instructive. We ignore the trivial case k, = 0. Since ^(O) = and k > 0, 
if K, is decreasing at any point y > 0, k,{z) must be decreasing V2; > y, since if 
not, concavity would obviously be violated. So suppose 3y 3 k{z) is decreas- 
ing V2; > y. Then we must have k{z) — ;► as 2; — > 00. Let yi > y. Then since 
K is decreasing on [y,oo) we must have K{yi) < K{y). Write rrii = HM^Jllyil^ Let 

1/2 > Vi — i^{yi)/f^ > Hi (recall m < 0). Then since K,{y2) > 0, 1712 := '^ "'^ < mi, 
and thus K,(ty + (1 — t)y2) > ti^iv) + (1 ~ 't)i^i3J2), so k is convex from y to y2- Thus 
we conclude that if /t(^) is concave with /t;(0) = 0, k must be nondecreasing. 



A. 3. 2 K concave ^ v continuous 

It suffices to prove that k is continuous, in light of the fact that the exponential 
function is continuous, and the composition of continuous functions is continuous. 
Suppose s < X < y < t. Now let Ai be such that Ais + (1 — Ai) y = x, and A2 be such 
that A2 a; + (1 — A2) t = y. Obviously as y — > x, Ai | 0, and A2 T 1- We have 

n{x) > Ai /t(s) + (1 — Ai)/t(i/) 
n{y) > A2/€(a;) + (l- A2)K(t). 

Since s and t are fixed, we have, as y ^ x, 

k{x) > lim {Ai/t(s) + (1 - Ai)k(i/)} = K{y) 

Alio 

k{x) < lim < T- K{y) + ^-^/t(t) \ = K{y). 

A2T1 I A2 A2 J 



Thus as 1/ ^ x, we must have n{y) -^ k{x). 
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A. 3. 3 v{x) = e^^l^l) is neither log-concave nor log-convex 

Certainly v restricted to certain subsets of its domain is log-concave, but V x, we have 

\og{v){h+{l-^){-x)) = K{0)=Q 
< k{\\x\\) = - \og{v){x) + (1 - 2) ^og{v){-x) 



violating log-concavity in general. 




Figure A.l: For admissible w, K(|a;|) is neither concave nor convex 



A. 4 The GRS-condition 

To illuminate further the GRS-condition, note that v{x) = e^'"' satisfies (3.2): 



l/n 



e — > 1 Va; G 



but Vo(x) = e'""' does not: 



r, \ l/n „ 



00 Vx G 



Note, however, that condition (3.1) does not imply condition (3.2). Let us restrict 
ourselves to d = 1 and treat only x > 0, to simplify the computations. Consider the 
example v = e'''^'^'^ where K(t) = t'"'^*-' and r is asymptotic to 1 from above: 



VAX 



3(X+1)^-1 
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We denote Ki(x) = (x + 1)^+t — 1. By construction f (0) = e'*'-*'^ = 1. The concavity 
of Ki can be shown by demonstrating that the second derivative k," is negative where 
x> 0\ 

But in checking the GRS-condition, we have 

(1 / \ -L / nx-\-2 \ 

ViinxV^'') = - (nx + 1)'^^^^^^^+^'' - 1 
^ n i 

1 1 1 

= —(nx + 1) ■ (nx + 1) ("^+1) 

n n 

1 

> X ■ {nx + l)("^+i) 

> X, (A.-24) 
since y^''^ > 1 if y > 1. Thus 

1 / \(nx+lV (nx+l)'_l] 

\imvi{nxY''' = hme'^ ^ ^ > e" > 1, (A.-24) 

n— +00 n— +00 

and so vi{x) does not satisfy the GRS-condition. 

A. 5 An example of a subconvolutive weight 

The admissible weight Ts{x) = (1+ jxl)'', s > 2d, is subconvolutive. To show this, fix 
X and divide M into two regions: 

t e K^ =^ It - xl < — 

I 1-2 

tern => |t-a;|>^. 
X 2 



For t G Nx, we have 



We conclude that 



ry /y /y -i /y 

<t-x<- => -<t< — 

2 - -2 2--2 

ltl>M ^ (i + |t|)-<(i + M)- 

:i + 1^1)"^ < (2 + 1^1). ^ 2^(1 + 1^1)"'- (A.-28) 



Thus 



{l + \t\)-%l + \x-t\y'dt < 2^(l + |x|)~' / {l + \x-t\ydt 



< 2^(l + |a;|)-* f {1 + \t\)-' dt. (A.-2^ 

Jr 



^I confess I have only shown this indirectly with Matheniatica; the inequalities involved are unsolvable in closed- 
form. The constrained maximum of k"{j/) with J/ > appears to be at ki (0.295403) = —0.419112. 
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On A^^, we have 1 + |t - a;| > 1 + ifi, so (1 + |t - x\)-' < (1 + ^-^y. Thus we have 
[ {1 + \t\)-%l + \x - t\)-' dt < (1 + M)- /■ (i + |t|)--rft 

< {1 + ^-^r f{i + \t\rdt. 

If s > 2, the integrals converge, and we have t^^ * t^^ < Ct^^ as was required. 

A. 6 Non-separability of Cl and Cy 

Neither C^ nor Cl is separable. To see this, consider the interval / = [0, 1]. Represent 
each z E I hj its binary equivalent .did2d^ . . . . Form the £°^(Z) sequence Hz via 
{yx)k = dk- "Unfold" Hz by defining z as follows: 

~ _ \i.yz)-2k+i j eZ, j < 0, 

[{yz)2k J ez, J >o. 

Then form a function fz = '^Xlk,k+i){x)xik,k+i){y)zk- Thus in each square centered 

kez 
diagonally on the line x = y, fz assumes only either of the constant values 1 or 0, 

and off the diagonal squares fz assumes only zero. 

Denote the collection of such functions by i3{o,i}- It is clear that i3{o,i} C Cy H Cl- 
We have \\fz^ - fzAc. > 1 and \\fz^ - fzAc\ > 1 V fz^,fz2 e >C{o,i} with zi ^ Z2 
since w(0) = 1 is a global minimum. So, in the topology generated by either norm, 
Bi{fzi) H-Bi (/^j) = 0. Since [0, 1] is uncountable, there exist an uncountable number 
of balls of radius |, containing only singletons, around the elements of -Cjo,!} (and 
£{0,1} is not even dense in £„ fl Cl). Thus £{o,i} is not separable, and hence neither 
is C,, nor Cl. 
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Appendix B 

Gabor analysis, LOFDM, and the 
Stone- von Neumann theorem 



B.l Powers of the Gabor frame operator commute with TF- 
shifts 

For all t G ]R\{1}, let T be the commutator 

T := S'T^M^ - T^M^S\ (B.O) 

Now S'^^*T5'^^* = 0, and since the frame operator S* is a positive operator, the 
operators S^ ^ u ^ 0, are well-defined and positive by the spectral theorem and thus 
invertible. We conclude that T = 0. In particular 5*^2 commutes with the time- 
frequency-shifts as well, and we have 

S-^g,A) = iS--^g,A). (B.O) 

B.2 Simultaneous pulse and lattice dilation commutes with 
Lowdin orthogonalization of the Gaussian 



We prove (7.1). First note that the Gram matrix of {g^, l/y^ap, y o j p) equals the 
Gram matrix for ((7, l/v^, 1/v^)- Consider an entry in the Gram matrix of dilated 
Gaussians on the dilated lattice: 

k9l.^,,^.. 9l^,,^) = /(2a)ie-'^-(*-^'^)^e-^'^^*(2a)ie---(-'=^)^-^^v^*rft. 
^ ^ Jr 

(B.O) 

By changing variables to r = ^/at, we get 
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From this we can conclude that for the Gram matrices derived from dilated and 
undilated Gaussians, we have 

leading to 

k,e k,e 

Now since the dilation operator "D is a bounded linear operator on L^(]R), we have 
from, e.g., Lemma 5.3.2 in [22] that 

Va = Yl ^ka%fifi 9ak,l = ^^ Yl ^ka/Pfifi 3k/ = V^ (/?. (B.O) 

B.3 Representations, The Stone- von Neumann theorem, and 
the Heisenberg group 

A representation ( oi a group G on a vector space V is defined as a homomorphism 

(■■G^GLiV), (B.O) 

with dim( = dimV, and where GL{V) denotes the group of linear operators on V. 
The rotation groups on M^ are examples of representations on a vector space without 
regard to a choice of basis. If we choose a basis, this defines an isomorphism 

GL{V) -^ GLniF), (B.l) 

T — > (matrix of T), (B.2) 

where GLn{V) denotes the group of invertible n x n matrices and F is the field over 
which the chosen basis lies. In this way we arrive at ( — > GL{V) — > GLn{F) 
which is a matrix representation of G. A matrix representation can be viewed as a 
representation of G on the space of column vectors F". 

Let C be a representation of G. Denote the image oi g E G under ( as (g. Then (g is 
a linear operator on V, and for g,h E G the associative law yields (gh = CgCh- 



Definition B.3.1 (Unitary Representations): A Hermitian vector space is defined as 
a complex vector space together with a positive-definite Hermitian form (■,■). Recall 
that an operator T on V is unitary <^ (x, y) = {Tx, Ty), Va;, y E V. It is thus natural 
to define unitary representations via {u,w) = {(guXgw), Wu,w EV,'ig E G. 



APPENDIX B - GABOR ANALYSIS, LOFDM, AND THE SVN THEOREM 95 

Consider the dihedral group D4 acting on M^; there are two proper invariant sub- 
spaces - the plane containing the square A on which D4 is defined, and the axis of 
rotation perpendicular to A. 

A proper G-invariant subspace of V under (^ is a subspace of V, not equal to V, which 
is invariant under the action of the representation (. 



Definition B.3.2 (Irreducible Representations): If a representation (^ of a group G 
has no proper G— invariant subspace, then ( is called irreducible. 



B.3.1 Schrodinger's Representation of the Heisenberg Group 

Define on H' = (R x M x T) a group multiplication by 

(x, y, r) ■ (m, v,ri) = {x + u,y + v, e2-(-+'?)e*-(«s'--)). (B.2) 

This is the reduced Heisenberg group. Further, let EI = (M x M x M) and define a 
multiplication 

(x, y, t)- {u,v,r]) = {x + u,y + v,T + r] + -{uy - xv)), (B.2) 

and the result is the full Heisenberg group. Both the full and reduced Heisenberg 
groups have the same distinguished unitary representation, the Schrodinger represen- 
tation. The Schrodinger representation vr from EI into the unitary operators on L^(IR) 
is given by 

7r(x,i/,r)/(t) =re"*^^e=^"^^V(t + a;), x,i/, t G M, r G T, (B.2) 

or 

7r(x,|/,r)=re-'^^M,T_,. (B.2) 

The Schrodinger representation is clearly unitary, and is irreducible since if we suppose 
/C 7^ is an invariant subspace of Ti under tt and fix (yf 7^ in /C and let / G /C"*", then 
the representation coefficients equal for all h = (x, y, r) G EI: 

0=|(/,7r(h)(?)| = |K,/(x,t/)|, V(x,|/,r)Ge (B.2) 

and since the STFT is injective, / = 0, so /C"*- = and )C = H. 

For convenience we will often write z = {x,y) for x,y E M."^ so n = 7i{z,t). Every 
lattice in M^ is of the form A = L{7?) for some L G G'L2(M). Since we use the 
tools of representation theory of the Heisenberg group, we must first determine which 
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automorphisms z h^ Lz oi M? extend to automorphisms ii{z,T) = {Lz,t) of H. If 
L leaves the symplectic form [{x,y), {u,v)] = yu — xv on M^ invariant, L defines an 
automorphism of EI in the following way: 



iiiz, t) = {Lz, r), V^ G M^ r G T, 
yielding the sequence of mappings 



e^e^^w^^:?^-^?^ 



Such a matrix L is characterized by^ 



det I L ■ 



y V 
jc u 

detL = ±1 

L G SL2{R) 



±det 



y V 
X u 



(B.2) 



(B.2) 



(B.3) 

(B.4) 
(B.5) 



Theorem B.3. 3 (Stone-von Neumann) Any irreducible unitary representation ( of 
EI on a Hilbert space Ti satisfying C(0, t)/ = r/, t ^T, f & 7i, is unitarily equivalent 
to the Schrodinger representation vr. Thus there exists unitary U : Ti. ^ L^(IR) such 
that C{x, y, r) = f/~%(z, r)f/, \/{z, r) G H. 

For L G S'L2(M) we compose the Schrodinger representation with the automorphism 
iL to yield a new representation ttl := n o i^ given by ttl^Zjt) = n{Lz,T). Clearly 
TiL is unitary, and is irreducible since L is a unitary change of coordinates, so (see 
the proof irreducibility for n above) the appeal to the STFT's injectivity still holds. 
Since 7rL(0,r) = 7r(0,r), Vr G T, then by Theorem B.3. 3, there exists unitary fi{L), 
called a symplectic operator, (defined only up to a phase factor) such that 



TiLiz, r) = fiiL)n{z, r)^(L)-^ y{z, r) G H. 



(B.5) 



The phase factors of the symplectic operators /i(-) can be defined so that either 
fj,{AB) = fi{A)fi{B) or fi{AB) = —fi{A)fi{B) (but neither case can hold universally), 
by invoking Mp{d), the two- fold cover of Sp{d). The following theorem is due to Shale 
and Weil ([50], [59]). 

Theorem B.3. 4 Let Mp{d) be the two-fold cover of Sp{d), and let A ^ A be the 

projection from Mp{d) onto Sp{d) = Mp{d)/{±I}. Then there exists a unitary rep- 
resentation fi of Mp{d) on L^(IR"') such that 



li{A)7i{z,T)fi{A) ^ = n{Az,T). 
We call Mp{d) the metaplectic group and /i the metaplectic representation. 

^In higher dimensions, this characterization defines the symplectic group. 



(B.5) 
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B.3.2 Examples relevant to Sections 7.2.3 and 7.2.4 

In this section, we drop the reference to r in (B.3.1), since it is of no consequence to 
the calculations and its inclusion would only be cumbersome. Note that each of the 
following matrices A, B, and C are members of 5*^(1). 



1 
-1 
(6.3.2), (6.3.2), and (B.3.1), we have 



Example B.3.5 Let A 

1 {B.c 



then IJ,{A) = T . Indeed, by equations (6.3.2), 






[B.5) 



Example B.3.6 Let B 



a 
a'^ 



; and set Vaf^x) = a^f^ f{ax). Then we have 



a 



a 



-l/2g2™(a-4-.)^(^-l^_^) 



-1/2 2m(a-'^Lo)(t-ax) 



f{a '^(t — ax)) 



T^.M^-i^V^-if{t). 



(B.4) 



Thus fi{B) =V^-i. 



Example B.3.7 Let C 

C^^ results from 



1 0' 

/3 1 



with (3 = (3, and set Cp = e'^'f^* f{t). Then ^i{C) 



C(}-^7r{x,uj)C^f{t) = e"*'^*'e"*^^e2"^'^(*-^)e-"*'^(*'^)V(i-a:) 

_ ^nix{uj+f3x) ^2ni{uj+f3x){t-x) r/^ _ ^\ 

= n{x,u; + Px)f{t) = n{C{x,uj))f{t). 



(B.3) 



B.3.3 Extension to general symplectic lattices 

The following is from [22]. In the context of frame theory, we have, for a window g, 
arbitrary / G -L^(]R), and an arbitrary lattice A 



Sg,Af = ^(/, 7r(A, T)g) 7r(A, T)g. 



(B.3) 



AeA 



Lemma B.3. 8 Suppose A = MI? ^ g G 5*0 (M) (functions whose STFT is integrahle) . 
Writing M = aL, det L = 1, Ai = (aZ)^ we have 



Sg,A = a{L)S^(^L)-^g,Ai'^i^) ' 



(B.3) 
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where a is the (unitary) metaplectic representation. 
Proof. Write S = Sg^A. Then 

Sif) = 5](/,7r(A)(?)7r(A)(? (B.4) 

AeA 

= 5^(/,7r(L«%)7r(L«% (B.5) 

= J2{f,ar{L)n{ak)a{L)-'g)a{L)n{ak)a{L)-'g (B.6) 

= a{L)J2HL)-'f,^ic,k)a{L)-'g)n{ak)a{L)-'g (B.7) 

fcez2 
= [a(L)5.(^)-.,,A,a(L)-i](/). (B.8) 



Thus the frame operator on a general lattice is unitarily related to the frame operator 
on a separable lattice. This fact permits us to conclude that all results on Gabor 
frames extend to the general lattice case. 
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Appendix C 

Glossary of notation 



In this appendix we provide a referral for the notation commonly used in this paper. 



C.l 


Basics 


n 


arbitrary Hilbert space 


Bin) 


the space of bounded linear operators on H 


n 


the nonnegative integers 


z 


the integers 


Zh. 


the positive integers 


Z'^ 


the set of rf-tuples of integers 


Q 


the rational numbers 


R 


the real numbers 


R+ 


the positive real numbers 


R'^ 


the set of li-tuples of real numbers 

V-1 

the complex numbers 


C 


y 


the complex conjugate oi y = a + hi: y = a ~ bi 


3 


"there exists" 


V 


"for all" 
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the empty set 

O the space of all functions on M^'^ 

C°° (M^) the space of all infinitely differentiable functions on M*^ 

Co {/ I / continuous and f{x) — ^ as |a;| -^ oo} 

\x\ the finite-dimensional 2-norm of x 

a (tti, . . . , ad) G Z'^; a multi-index 

x" nil X? 



op 



\ dxi I \ dxi 



S the Schwartz space; all ip G C°^(]R'^) such that 

sup \x°' D f^ (p{x)\ < oo 



S* the space of tempered distributions, the dual of S 

5 the Dirac delta function (point measure) 

!F the Fourier transform operator 

STFT the short-time Fourier transform 

Af the ambiguity function of / 

A{f, h) the cross-ambiguity function of / and h 

f -k g composition of integral operators with kernels / and g, also 

frequently denoted fg 

f * g the convolution of / with g 

W{W^) the Wiener algebra : / G W{W^) ^ /, / G L^i.W^) 

f f{x) dx the Lebesgue integral of the function / 

the norm of / 

the operator norm of / given by sup ||/(a;)|| 



llx!l=i 



the p-norm of /, given by ||/||p = {J^^ \f{x)\Pdx)'' 
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L (W^) the space of p*^-power Lebesgue-integrable functions on M*^, 

1 < p < oo 

ess sup f{x) the essential supremum of /, given by 

ess sup f{x) = inf {M | m{x \ f{x) > M} = 0} 



the infinity-norm of /, given by ||/||oo := ess sup |/(a;)| 

\f the Frobenius norm of the matrix A, given by Yl^i Yll=i k«iP 

A"^ the complement of A in some superspace 



C.2 Banach algebras and off-diagonal decay 

S{x — y) the kernel of the identity integral operator on / : M'^ — > M; 

also denoted e 

Q the space of all functions on M^*^ appended with 6{x — y) 

0,1 a superspace of Q admitting point measures, of which no 

hyperplane can contain an infinite number 

V an admissible weight, which satisfies the GRS-condition and 

whose exponent n satisfies concavity on |x| 

K the function in the exponent off : v{x) = e''*-'^'^ 

f Q g the Hadamard (or pointwise) multiplication of / and g 

jCI the space of equivalence classes of all operators in Qi which are 

integrable against the weight v in both coordinates 

jC^ Cl with v = l 

Cy the space of equivalence classes of all operators in fii which are 

ess-sup-bounded pointwise against v in both coordinates 

Ts the polynomial weight (1 + |a;|)* with s > 

Bu,s derived from C^ : v = UTg with u admissible; Bu,s = /^^ fl £„ 

/i the convolution (Laurent) operator generated from / : M^ ^ M 

/ the non-distributional part of the kernel / of an integral operator 
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PaU) 

pU) 



the algebraic spectrum of the operator / 

the spectrum (as an operator on Ti) of the operator / 

the algebraic spectral radius of the operator / 

the spectral radius (as an operator on Ti) of the operator / 



C.3 PseudodifFerential operators 



H 

a 

A 



time-varying communications channel operator 

a pseudodifFerential operator 

the symbol of K^ 

the kernel of K„ 

the partial Fourier transform in the second variable 

a linear time-varying system 



C.4 OFDM and LOFDM 



OFDM orthogonal frequency- division multiplexing 

BFDM biorthogonal frequency-division multiplexing 

EBFDM extended biorthogonal frequency- division multiplexing 

SNIR signal-to-noise-plus-interference ratio 

AWGN additive white Gaussian noise 

TFL time-frequency localization 

T the symbol length for a rectangular OFDM system 

F the carrier separation for a rectangular OFDM system 

p the spectral efficiency of a communications system 

Ch scattering function of the channel H 

S'h spreading function of H (randomized version of Ch) 
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h{t, t) time-varying impulse response of H 

Tq maximum delay due to multipath propagation 

uq maximum Doppler spread 

At/ second-order time-localization of the function / 

Ai/f second-order frequency- localization of the function / 

■0 the atom, or generating pulse-shape, of an OFDM system 

/fe^^ f{t — A;T)e^'^*^-^*; notation used only for rectangular OFDM systems 

Tx the translation operator: Txf{t) = f{t — x) 

M^ the modulation operator: M^f{t) = f(t)e^'"^'^-* 

z{\) T^M^:\ = {x,iu)eA 

A an arbitrary symplectic lattice 

L the generator matrix for the lattice A 

^(A) the density of the lattice A 

Ao the adjoint lattice of A given by ^^^^^A 

J an arbitrary index set 

Cj a frame element 

S the general frame operator: Sf = J2jej{fy^j)^j 

S'^{g, A) the action of the z/*^-power of the Gabor frame operator 

S on the entire Gabor frame {g, A) 

Sf^h,A the frame operator for the frame generated by the dual functions 

/ and h on the lattice A 

(/, A) Gabor system generated by / on the lattice A; also used to 

denote an LOFDM system 

Lo(/, A) The Lowdin-orthogonalized pulse derived from / on an 

arbitrary time-frequency lattice A 

(/, T, F) same as (/, A) but with A rectangular with parameters T and F 

Ld{f, T, F) Ld{f, A) where A is rectangular with parameters T and F 

Ld{f,A) 
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</7o- Ld{gfj, A) where g„ is a dilated Gaussian 

9? Ld{g„,A) where a = 1 

supp^f the effective support of /, where |/| > e 

Rf^A the Gram matrix of (/, A) 

Rf,T,F the Gram matrix of (/, T, F) 

Vafit) the dilation operator ^/af{at), a G IR+ 

Cfif{t) multiphcation by a chirp: f{t)e-'^'^^\ /3 G M 

5*^(2, M) the special linear group; 2x2 real matrices A with |det(A)| = 1 

C.5 Future research and appendices 

W a weight matrix under the usual matrix multiplication 

W° a weight matrix under Hadamard matrix multiplication 

B^{f) the ball of radius e about a function / 

[a, h] the set of real numbers x such that a < x < b 

V a vector space 

C a reperesentation of a group G on K 

^g the image oi g E G in the vector space V under the representation ( 

T the torus: {{x,y) eM."^ \0 < x,y < 1} with {0,y) and (1,1/) identified 
and (x, 0) and (a;, 1) identified 

W (M X M X T) 

e (M X R X M) 

TT the Schrodinger representation TT{x,y,T) = re^^^^MyT^x, 

x, y, t G M, r G T from either M.' or EI to the unitary operators 
on L2(M'^) 

Tiiiz, r) the representation 7r(Lz, r) V (2;, r) G EI for L G S'L2(IR) 

Sp{d, M) the symplectic group, also denoted Sp{d) 

Mp{d, M) the metaplectic group, the twofold cover of Sp{d) 
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//(•) the symplectic (unitary) operator appearing in the metaplectic 

representation of Mp{d, M) 
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